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ABSTRACT 

This  paper  considers  a  joint  replenishment  inventory  problem  with  a  continuous- 
review  (S,  c,  s)  policy  for  the  backorder  case  with  Poisson  demands  and  constant 
procurement  lead  times. 

Whenever  item  i's  inventory  level  hits  sj  (reorder  point)  or  lower  it  triggers  an 
order  so  as  to  raise  item  i's  level  to  S-  (order  up  point).  At  the  same  time  any  other 
item  j  with  inventory  level  at-or-below  its  can-order  point  C:  is  included  in  the 
replenishment. 

A  Poisson  demand  model  with  a  queueing  description  of  the  system's  operation  is 
analysed,  and  comparisions  are  conducted  for  joint  versus  individual  orders  in  the  case 
of  multi-item  problems,  where  joint  replenishment  of  several  items  may  reduce -setup 
costs. 
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I.  INTRODUCTION 

A.       DESCRIPTION  OF  THE  PROBLEM 

The  purpose  in  constructing  a  mathematical  model  of  an  inventory  system  is  to 
use  it  as  an  aid  in  developing  a  suitable  operating  doctrine  for  the  system.  The 
criterion  most  frequently  used  for  selecting  the  operating  doctrine  is  that  of  profit 
maximization  or  cost  minimization.  In  some  cases  the  task,  of  determining  the  optimal 
operating  doctrine  is  so  difficult  that  it  is  either  impossible  or  uneconomical  to 
determine  the  optimal  doctrine,  and  instead,  one  optimizes  with  respect  to  some  subset 
of  operating  doctrines.  Occasionally,  the  mathematical  model  may  be  so  complicated 
that  it  is  extremely  difficult  to  do  anything  analytically.  In  such  situations,  simulations 
are  used  to  study  various  operating  doctrines. 

We  are  concerned  with  finding  optimal  ordering  policies  for  a  multi-item 
inventory  system  in  which  the  inventory  position  is  under  continuous  review.  In  this 
thesis,  the  times  at  which  demands  occur  are  assumed  to  be  generated  by  independent 
Poisson  processes  with  intensities  X-  (annual  demand  rate  for  item  i).  We  assume  that 
all  shortages  are  backordered  and  will  be  satisfied  from  the  next  shipment. 

We  assume  that  a  holding  cost  of  H-  per  unit  time  is  charged  for  each  unit  of 
item  i  in  inventory.  When  an  order  is  placed,  a  setup  cost  for  an  independent 
replenishment  of  item  i  is  assumed  to  be  C-  =  A  +  aj,  where  A  is  a  fixed  cost  and  a-  is 
a  cost  which  depends  on  the  item  ordered.  The  cost  of  a  joint  replenishment  of  both 
item  i  and  j  is  A  +  a-  +  a:  .  For  instance,  if  item  i  triggers  a  replenishment,  the  setup 
cost  of  A  +  a-  is  applied  to  item  i,  and  a:  to  the  jointly  replenished  item  j  . 

This  type  of  cost  structure  is  particularly  appropriate  when  a  group  of  items  is 
ordered  from  the  same  supplier,  or  if  a  group  of  items  uses  the  same  means  of 
transportation,  and  transportation  costs  do  not  increase  proportionally  when  the 
quantity  transported  is  increased.  Thus,  when  a  replenishment  is  made,  there  is  a 
major  fixed  cost  independent  of  the  number  of  items  involved  and  a  variable  cost  that 
depends  upon  the  number  of  items  [Ref.  1:  pp.  278-83]. 

Because  of  this  cost  structure,  it  is  clear  that  there  may  be  an  opportunity  to 
reduce  total  variable  costs  if  one  takes  advantage  of  joint  replenishment. 


B.       SCOPE 

Our  interest  in  this  problem  stems  from  the  "random  joint  order  policy''  proposed 
by  Balintfy  [Ref.  2],  and  from  an  (S,c,s)  policy  treated  by  Silver  [Ref.  3].  For  the 
backorder  case  with  Poisson  demands  and  constant  lead  times,  and  the  specified 
customer  service  levels,  we  will  compare  the  individual  order  policy  versus  joint  order 
policy. 

However,  for  selecting  the  control  parameters  of  (S.c,s)  in  Joint  Order  Policy,  we 
have  considered  slightly  different  procedures  suggested  by  Schacck  and  Silver  [Ref.  4], 
and  Silver  [Ref.  3]. 

Two  programs  were  written  in  order  to  conduct  the  simulation  experiments.  The 
first  program  determines  the  control  parameters  of  each  item  in  the  group  for  both 
independent  and  joint  replenishment  policies.  A  second  program  evaluates  the  model. 
Detailed  user  instructions  and  the  program  itself  are  contained  in  Appendix  A  and  B, 
while  examples  of  their  use  are  presented  in  Chapter  II. 

The  individual  order  policy  is  demonstrated  only  for  the  purpose  of  comparison. 
For  the  sake  of  simplicity,  we  compare  the  joint  order  policy  to  the  individual  order 
policy  where  the  classical  lot  size  formula  is  used  for  each  item. 


II.  GENERAL  FORMULATION  OF  THE  MODEL 

A.       ASSUMPTIONS  AND  NOTATIONS 

1.  Assumptions 

We  are  dealing  with  a  group  of  items  where  the  cost  of  replenishing  two  or 
more  items  at  the  same  time  is  less  than  the  total  cost  of  replenishing  the  same  number 
in  separate  individual  replenishments.   The  assumptions  of  the  model  are: 

1.  There  is  fixed  cost  A  associated  with  each  replenishment,  and  variable  cost  a, 
associated  with  each  item  involved  in  the  replenishment. 

2.  Demands  for  item  i  are  Poisson  distributed  with  parameter  Xj. 

3.  There  are  no  quantity  discounts. 

4.  Inventorv  holding  costs  are  proportional  to  the  average  dollar  value  of 
inventory. 

5.  Service  level  is  defined  in  one  of  two  wavs:  probabilitv  of  no  shortase  per 
replenishment  cvcle,  and  fraction  of  demand  to  be  satisfied  directlv  from  the 
shelf. 

2.  Basic  Notation 

We  use  the  following  notation: 
A   :  fixed  cost  per  replenishment,   (independent  of  the  numbei  of  items  and 

units  involved) 
a-  :  fixed  cost  for  including  item  i  in  a  replenishment. 
Cj :  ordering  cost  per  order:  (  Cj  =  A  +  a-) 
EC-:  expected  relevant  annual  inventory  cost  for  item  i. 
Pj    :  desired  customer  service  level,  measured  in  terms  of  the  probability  of  no 

shortage  per  replenishment  cycle. 
P2   :  desired  customer  service  level,  measured  in  terms  of  the  fraction  of 

demand  satisfied  directly  from  stock  on  hand. 
L    :  replenishment  lead  time  in  years, 
r    :  inventory  holding  cost  rate, 
v^  :  standard  unit  price  of  item  i. 
s^ :  reorder  or  must-order  point  of  item  i. 
c-  :  can-order  point  of  item  i. 
S-  :  order-up-to  point  of  item  i. 
Qj  :  order  quantity  of  item  i. 
X-  :  Poisson  demand  rate  for  item  i,  in  pieces  per  year. 
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I:   :  average  on  hand  inventor/  level  of  item  i. 
n      :  number  cf  items  in  the  replenishment  group. 
N:   :  expected  number  of  replenishments  which  involve  item  i  per  year. 
NT:  :  expected  number  of  replenishments  triggered  by  item  i  per  year. 
Other  notation  will  be  introduced  as  they  needed. 

B.       INVENTORY  THEORY  BACKGROUND 

Before  considering  the  coordinated  control  problem  it  might  be  helpful  to  briefly 
discuss   the   related   single   item  models  whose   solution  will   be   a   key   element   in 
determing  the  parameters  of  the  (S,  c,  s)  control  system. 
1.  Economic  Order  Quantity  Model 

The  E.O.Q.  model  is  based  on  the  assumption  that  the  entire  lot  is  added  to 
stock  at  one  time,  and  that  the  stock  will  be  withdrawn  at  a  constant  rate,  and  no 
stockouts  are  permitted. 

Total  annual  variable  cost  =  Ordering  cost  +  Holding  cost 

EC£-  Xiq/.Qi+  HjQr/2  (2.1) 

The  annual  order  cost  is  obtained  by  multiplying  the  number  of  orders  per 
year  (X.-/Q-)  by  the  cost  of  placing  an  order  (C-). 

The  average  annual  holding  cost  for  item  i  is  the  average  inventory  (Q:/2) 
times  the  annual  unit  holding  cost  (H-=  v-r). 

To  obtain  the  minimum  cost  lot  size  (E.O.Q. ),  take  the  first  derivative  of  total 
annual  cost  with  respect  to  Q-  and  set  it  equal  to  zero.  Solving  this  for  Q-,  we  get  the 
E.O.Q.  formula: 

Qi*  =  V  [2CiXi  /  Hj]  (2.2) 

Once  the  economic  order  quantity  is  known,  the  expected  number  of  orders 
placed  during  the  year,  n.  can  be  determined: 
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2.  Exact  Model  for  the  Backorders  case  with  Poisson  Demands  and  Constant  Lead 

Times 

If  demands  are  random,  the  probability  of  stockouts  exists,  and  we  must 
include  a  stockout  cost  in  the  expression  for  expected  annual  variable  costs.  This 
expression  is  given  in  [Ref.  5:  eqn.  (4-61)]  to  be 

EC-  *iCi/Qi+  H{  (Qj  /  2  +  1/2  +  sj  -  M()  +  ^  E(Qif  sj)  (2.3) 

+  (tj  +  H^Qi,  sj) 

where 

(p-  :  the  cost  of  backorder  per  unit  quantity  of  backorder. 

t-  :  the  cost  of  backorder  which  is  proportional  to  the  length  of  backorder  time. 

M:  :  expected  lead  time  demand  for  item  i  in  pieces,  M-=  X-L. 

H-  :  holding  cost  of  item  i. 

E(Q-,  S:)  :  expected  number  of  backorders  incurred  per  year  for  item  i. 

B(Qj,  s[)  :  trie  steady  state  expected  number  of  backorders  for  item  i. 

It  is  not  easy  to  derive  the  optimal  Q-  and  s-  by  use  of  the  above  exact 
expression  for  EC-.  Fortunately,  it  turns  out  in  practice  that  it  is  seldom  necessary-  to 
use  the  exact  formulation,  except  for  the  case  when  it  costs  very  little  to  incur 
backorders.  The  approximation  discussed  in  the  next  section  has  been  shown  to  give 
good  results  when  backorders  are  infrequent. 

3.  Approximate  Form  with  Known  Backorder  Cost 

The  average  annual  variable  cost  equation  can  be  approximated  as  follows 
[Ref.  6:  pp.  59  -  62]: 

ECr  Xiq/Qi.+  Hi(Qj/2.+  si-  M^  +  cp^i/Q^MpSi)  (2.4) 

where 
q>-  :  the  cost  of  backorder  per  unit  quantity  of  backorder. 
E(M->  s-)  :the  expected  number  of  backorders  at  the  end  of  cycle. 

To  minimize  EC-,  take  the  partial  derivative  of  the  ECj  with  respect  to  Q-  and 
s-  and  set  these  equal  to  zero. 

dEq/dQ-  -  ^q/Qi2  +  Hj/2  +  q»iX.iE(Nfi>si)/Qi2-  0  (2.5) 
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dECj/dsj-  Hi  -r  (f-Xj  /  Q^dECM^  s^/ds-  0  (2.6) 

The  first  equation(2.5)  yields  the  optimal  Qj  for  a  given  reorder  point  s-: 
Qi*  =  V  [2Xi(Ci  +    (PiECM^  sj)  /  HJ  (2.7) 

The  second  equation(2.6)  yields  : 

dECM^/d^  =  -JP^Mi)  (2.8) 

where 

E  (Mj>  s-)   :  the  expected  number  of  backorders  at  the  end  of  cycle. 
-  Y,  Pr  (Mi)  :  complementary  cumulative  distribution. 

This  is  the  complementary  cumulative  distribution  of  the  Poisson  random 
variable,  X,  evaluated  at  s-.  Therefore,  solving  the  second  equation  for  s-  in  terms  of 
Q-  gives  the  maximum  allowable  probability  of  a  shortage  during  a  lead  time  : 

F'(Si)  =  (H^yC^i)  (2.9) 

To  find  the  optimal  pair  (Q^,  s^)  that  minimize  EC^,  iterative  procedures  can  be 
applied,  such  as  described  in  [Ref.  6:  p.  61].  Hower,  we  do  not  apply  this  kind  of 
procedure  to  find  out  the  values  of  control  parameters  in  (S,  c,  s)  policy,  because  the 
order  size  in  the  (S,  c,  s)  policy  is  not  fixed. 

In  (S,  c,  s)  policy,  the  order  size  is  determined  as  an  order-up-to  point  minus  a 

reorder  point  for  the  item  which  triggers  the  replenishment  action.    Order  sizes  for 

jointly  replenishable  items  are  order-up-to  points  minus  on-hand  levels  which  are  less 

than  or  equal  to  can-buy  points  but  more  than  reorder  points. 

4.  Some  Useful  Formulas  for  Obtaining  the  Imputed  Backorder  Cost  /  Reorder 
Point  /Safety  Quantity  in  Fixed  Order  Size  System 

In  any  practical  situation,  it  is  very  difficult  to  determine  accurately  backorder 

costs.    They  can  include  such  factors  as  loss  of  customers'  goodwill(i.e.,  in  the  future, 

he  may  take  his  business  elsewhere),  or  in  the  military  supply  system,  the  cost  of 

having  some  first  line  weapon  system  inoperative  because  of  lack  of  parts.     Other 

contributions  to  shortage  cost  can  be  somewhat  easier  to  measure.     However,  these 

are  usually  small  part  of  the  total  backorder  costs  [Ref.  5:  p.  18]. 
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We  shall  assume  here  that  there  is  a  backorder  cost  associated  with  each  unit 
backordered. 

A  complete  cost  minimization  analysis  of  the  optimal  values  of  the  model 
parameters  can  be  made  only  if  a  cost  can  be  assigned  to  each  shortage.  As  an 
alternative  to  specifying  the  shortage  cost,  one  can  examine  a  given  inventory  policy 
and  obtain  the  imputed  cost  of  shortage  within  it  [Ref.  7:  p.  338].  Rather  than 
explicitly  costing  backorders,  we  will  specify  a  service  level  constraint  and  determine 
the  imputed  stockout  cost. 

Whenever  an  organization  uses  a  service  level  to  establish  a  reorder  point,  it 
effectually  establishs  a  stockout  cost.  Associated  with  a  given  service  level  is  an 
imputed  or  implicit  stockout  cost.  It  is  a  simple  matter  to  determine  the  imputed 
stockout  cost  for  a  given  service  level  from  the  probability  of  a  stockout  [Ref.  8:  p. 
162]. 

a.   Probability  of  no  Shortages  per  Replenishment  Cycle 

A  service  level  based  on  the  frequency  of  service  per  order  interval  or 

replenishment  cycle  will  indicate  the  probability  of  not  running  out  of  stock  during  the 

replenishment.    This  approach  does  not  concern  itself  with  how  large  the  shortage  is, 

but  only  with  how  often  it  can  occur  during  the  lead  time  (replenishment  cycle).    It  is 

defined  as  the  fraction  of  the  replenishment  cycles  without  depletion  of  stock: 

number  of  order  periods  with  a  stockout 

Service  level  fraction  per  cycle  =  1  -  

total  number  of  order  periods 
=  l-P/M^Si) 

Pj-CM^Sj)  =  Pr(Si) 

=   1  -  (service  level  fraction  per  cycle) 

The  term  Pf(s-)  is  the  stockout  level  fraction  per  order  cycle  or  the 
probability  of  at  least  one  stockout  while  awaiting  a  supplier's  delivery.  It  is  also  a 
measure  of  the  fraction  of  lead  time  periods  during  which  the  demand  will  exceed  the 
reorder  point.   The  magnitude  of  the  stockout  is  ignored  with  this  approach. 

The  following  formula  is  developed  to  compute  the  imputed  backorder  cost: 

Pr(Si)  =  (HjQj)  /  (bfa)  (2.10) 

where 

Pr(Sj):  probability  of  stockout  for  item  i. 
H-  :  holding  cost  of  item  i,  (H-=  v^r). 
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b;  :  imputed  backorder  cost  of  item  i. 

l 

To  obtain  the  appropriate  stockout  cost,  solve  eqn.(2.10)  for  b:.  An 
example  will  illustrate  the  procedure  in  a  later  section. 

b.   Fraction  of  Demands  Satisfied  Directly  from  the  Shelf 

The  fraction  of  units  demanded  and  immediately  filled  from  the  shelf  can  be 
defined  as  the  ratio:  (number  of  units  supplied)/(total  number  of  units  demanded) 
which  is  equal  to  (1  -  [(number  of  units  short)/(totai  number  of  units  demanded)]}. 

The  above  relationships  must  be  measured  over  some  fixed  time  period 
such  as  the  order  interval. 

For  example,  to  obtain  the  stockout  level  fraction  for  units  demanded: 


Stockout  level  fraction  for  units    = 


E(Mi>si) 


quantity  demanded  during  a  cycle 


E(Mi>si)/QL=  (HiQj)  /  (bjXj)  (2.11) 

where 
E(Mj>  s-):  expected  number  of  backorders  for  item  i  during  the  cycle. 
H-  :  holding  cost  of  item  i,  (H-  =  v-r). 
bj :  imputed  backorder  cost  for  item  i. 
c.  Safety  Stock 

With  backorders,  there  is  no  loss  of  sales  since  the  customer  waits  for  the 
arrival  of  the  next  order  if  stock  is  not  available.  The  expected  safety  stock  is  defined 
as: 


Safety  stock  =  £  (s  -  M)  Pf  (M)  (2.12) 

=  s  y  Pr  (M)  -  V  M  Pr  (M)  =  s  -  E(M) 
where 

E(M)  :  the  expected  lead  time  demand. 

The  number  of  backorders  per  lead  time  =  0,  if  M  <  s. 

=  M  -  s,  if  VI  >  s. 
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So,  the  safety  stock  is  simply  the  reorder  point  minus  the  average  lead  time 
demand.  The  reorder  point  is  determined  so  that  the  customer  service  objective  is 
satisfied. 


C.       MATHEMATICAL   STATEMENT   OF   THE    DECISION    PROBLEM    FOR 
THE  JOINT  REPLENISHMENT  MODEL 

1.  Cost  Equation 

The  basic  cost  equation  in  the  joint  replenishment  model  is  as  follow  [Ref.  3]: 

Eq=  I^r   +  NTjA  +  Njaj  (2.13) 

The  first  term  represents  the  inventory  holding  costs  and  the  other  two 
represent  the  set  up  charges  allocated  to  item  i.  The  total  expected  costs  per  year  of  the 
group  of  items  are  given  by 

EC  =  T  Eq=  Y  (Ijvjr    +  NT^A  +  Nja£)  (2.14) 

2.  Service  Level  Constraint 

The  service  constraint  takes  one  of  two  forms  depending  upon  the  measure  of  service 
used. 

(1)  Probability  of  no  shortage  per  replenishment  cycle 

This  is  simply  the  probability  that  demand  for  item  i  during  the  lead  time  will 
exceed  the  reorder  point. 


Pr(s)  =  Pr(M1>Sl)  =  EF-rflVii)  (2-15) 

where 

Pr(s):  probability  of  stock  out. 
M:     :  lead  time  demand  for  item  i. 

(2)  Fraction  of  demand  satisfied  directly  from  the  shelf 

The  expected  stockout  quantity  during  the  lead  time  is  given  by 

E(M  >  Si)  =  EdMi-^PrCMi)  (2.16) 
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where 

E(M:>  S:)  :  expected  number  of  stock  out  for  item  i  in  pieces 
during  the  lead  time 


*e 


M^  :  lead  time  demand  for  item  i 

A  backorder  occurs  if  and  only  if  the  lead  time  demand  (which  is  Poisson  with 
parameter  M-)  is  greater  than  the  inventory  level  of  the  item  when  the  order  is  placed. 

D.       INDIVIDUAL  ORDER  POLICY 

Under  the  the  individual  order  policy,  there  is  a  fixed  order  quantity  for  item  i 
that  is  ordered  every  time  the  reorder  point  is  reached  . 

Safety  stock  is  needed  to  protect  against  stockouts  between  the  time  the  reorder 
point  is  reached  and  the  order  received.  As  presented  in  [Ref.  3],  the  individual  order 
policy  (called  independent  control)  is  a  special  case  of  the  joint  replenishment  policy  (S, 
c,  s)  where  c-=  s-.  In  this  case  item  i  is  ordered  only  when  the  inventory  position  hits 
its  reorder  point,  and  the  order  quantity  is  always  of  size  Q-  =  S-  -  s-.  No  other  items 
will  be  included  in  the  replenishment,  hence  a  set  up  cost  of  A  +  a-  is  incurred  with 
each  replenishment. 

As  mentioned  in  the  introduction,  the  main  focus  of  this  thesis  is  to  compare 
joint  replenishment  with  independent  control  to  determine  the  cost  savings  that  can  be 
achieved. 

Let  us  first  consider  independent  control.  Then,  c-  =  s-,  and  the  optimization 
problem  is: 

Minimize  Eq  =  (Sj  -  Sj  +  1)/  2  +  (sj  -  J^Ltyr  (2.17) 

+  (A  +  a^  /  (S[  -  Si) 
subject  to 

Pr(X  ^  ti\X-L)  >  Pj  (2.18) 

or 

1-E(X>  s^^L)  'Q;>  P2  (2.19) 
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where,  X  is  the  lead  time  demand  with  a  Poissen  distribution  having  mean  XjL. 

The  left  hand  side  of  eqn.(2.1S):  Pr(X  ^  s-|  X^L)  is  the  probability  of  no  shortage 
per  replenishment  cycle. 

In  the  left  hand  side  of  eqn.(2.19),  the  term  E(X>  S:)/Q:  is  the  stockout  fraction 
of  units  demanded.  Therefore  1  -  E(X>s-)/Q-  is  the  fraction  of  demand  satisfied 
directly  from  the  shelf. 

The  exact  solution  to  this  problem  requires  that  the  pair  (Sj,  s-)  which  minimize 
eqn.(2.17)  subject  to  either  eqn.(2.18)  or  eqn.(2.19)  be  found  for  each  item  i.  However, 
it  have  been  observed  by  many  researchers  (see:  for  example  Silver  and  Wilson  [Ref.  9]) 
that  a  much  simplier  sequential  approach  does  almost  as  well  from  the  standpoint  of 
minimizing  costs.   This  sequential  approximation  approach  is  what  we  describe. 

1.  Applied  Solution  Procedure 

Stepl  :  Let  Q-  =  Sj  -  s-.  Find  the  value  of  Q-  which  minimizes  EC-  ignoring 
the  service  level  constraint.  Then  given  this  value  of  S-  -  s-,  the  constraint  is  used  to 
find  the  lowest  allowable  value  of  s-.  This  first  step  leads  to  the  usual  EOQ  expression 
Qi=  Vsi=  V[2(A  +  ai)^  HJ. 

Step2  :  Determine  the  reorder  points. 

(l)For  the  service  measure,  probability  of  no  shortage  per  replenishment  cycle, 
find  the  smallest  value  of  s-,  which  satisfies  the  condition  of  eqn.(2.18): 

Pr(X  <;  SilXjL)  >  Pj 

(2)For  the  service  measure,  fraction  of  demands  satisfied  directly  from  the 
shelf,  find  the  smallest  value  q!l  sj,  which  satisfies  the  following  inequality: 

XjL  -  Si  -  X{L  Pr  (X  <  Si  -  1  \\L)  +  SiPr  (X  <  s{  |  XjL)  ^  Q{{\  -  P2) 

Step3  :  Set  Si  =  Sj  +  Q-v 

Step4  :  Compute  the  average  annual  variable  costs. 

EC^  (S|  -  sj  +  1)/  2  +  (sj  -  X^vjr  +  (A  +  aj)3tj  /  (Sj  -  sj)  using  the  values 
Sj  and  S:  estimated  from  Steps  2  and  3.  A  computer  program  for  solving  for  these 
parameters  is  provided  in  Appendix  A. 

Note  that  we  can  now  obtain  imputed  stockout  costs. 
Step5  :  Compute  the  imputed  stockout  cost 
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(l)For  the  service  measure,  probability  of  no  shortage  per  replenishment  cycle: 
Use  eqn.(2.10):  Pf(s)  =  (H^)  /  (bjXj)  and  solve  for  b{,  bj-  HjQj  /  (1-  Pj)  Xy 

(2)For  the  service  measure,  fraction  of  demands  satisfied  directly  from  the 
shelf,  it  is  necessary  to  solve  for  E(M  >  Sj)  and  find  the  associated  Pf(s)  before  using 
eqn.(2.11).  It  can  be  done  alternatively  by  finding  the  complementary  cumulative 
distribution  of  the  lead  time  demand,  introduced  by  eqn.  (2.8).  Then  eqn.(2.11)  yields 
IT  H^M-Vp^Mi)]/^. 

2.  Numerical  Examples 

Repeated  from  Silver's  example  in  his  paper  [Ref.  3]. 


cvcle. 


Example   1.    Service  Level   =    Probability  of  no  shortage  per  replenishment 

Consider  the  following  example  involving  four  items: 
A  =  S  50 

a;  =  S  10  fori  =  1,  2,  3.4 
r  =  0.2  S  /  S  /  yr. 
L  =  1  month  (i.e.  1/12  year  ) 

TABLE  1 
INPUT  DATA 


item 
(i) 

(piecesvyr) 

V; 

(S;  piece) 

P, 
(desired  service  level) 

1 
2 

3 
4 

290 
41 

77 
122 

6.90 
1.20 
3.90 

2.30 

0.95 
0.95  . 
0.95 
0.95 

Use  of  eqn.  (2.15),  (2.17),  (2.18)  gives  the  results  shown  in  Table  2. 

Example  2.   Service  Level  =  Fraction  of  demand  satisfied  directly  from  the  shelf. 

Consider  the  same  data  as  example  1,  but  with  a  new  the  measure  of  service 
level.  Suppose  the  specified  fraction  of  demands  to  be  satisfied  directly  is  0.99.  The 
results  shown  in  Table  3. 
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TABLE  2 
INDEPENDENT  ORDER  POLICY:  SERVICE  LEVEL  95%  (Pj) 


Item 

LT  demand 

EOQ  (integer) 

Reorder 

TVC(S) 

SO  cost 

1 

2 

3 
4 

24.167 
3.417 
6.417 

10.167 

158.800(159) 
143.17SU43 
108.840(109) 
178.399  (178) 

33 

7 

11 

16 

232.024 

i;    1,11 

88.860 

S4.977 

15.113 

16.762 
22.051 

13.453 

TOTAL : 

441.203 

TABLE  3 
INDEPENDENT  ORDER  POLICY:  SERVICE  LEVEL  99%  (P->) 


Item 

LT  demand 

EOQ  (integer) 

Reorder 

TVC(S) 

SO  cost 

1 
2 
3 
4 

24.167 
3.417 
6.417 

10.167 

158.800(159) 

143.178(143) 
108.840(109) 
17S.399    178) 

25 
3 
7 

10 

220.984 

34.383 
85.740 

82.217 

"1.983 

1.SS2 

3.499 
1.536 

TOTAL  :     423.323 

E.       JOINT  ORDER   POLICY 
1.  Introduction 

Let  us  now  consider  the  possibility  of  joint  replenishments.  As  discussed 
earlier,  we  can  reduce  ordering  setup  costs  by  combining  the  orders  of  several  items 
whenever  an  order  is  placed.  This  may  result  in  increased  holding  costs  but  the  savings 
in  ordering  costs  may  be  greater  than  the  increase  in  holding  costs  resulting  in  a  net 
decrease  in  total  costs.  Further,  it  is  likely  that  the  occurance  of  stockouts  will  be 
reduced  since  we  will  frequently  reorders  some  items  before  the  inventory  levels  hit  the 
normal  reorder  points  for  those  items. 

We  do  not  seek  to  determine  the  "optimal"  joint  replenishment  policy. 
Rather,  we  restrict  attention  to  a  simple,  intuitively  appealing,  policy  and  seek  to  find 
the  optimal  values  of  parameters  for  the  selected  policy.  The  policy  we  consider  is  the 
(S,  c,  s),  "can  buy,  must  buy"  policy  which  requires  an  order  to  be  placed  for  item  i  up 
to  the  level  Sj  any  time  the  inventory  position  for  item  i  reaches  or  falls  below  s-.  In 
addition,  however,  we  order  any  other  item  i  whose  inventory  position  is  below  the  can 
buy  point  c-,  we  always  order  up  to  the  level  S-. 
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2.  General  Concepts  for  Obtaining  the  Values  of  Control  Parameters 

Suppose  a  single  item  i,  is  in  a  group  of  coordinated  items  and  assume  initially 
that  the  replenishment  lead  time  is  zero.  The  assumption  of  unit  sized  Poisson 
demands  together  with  a  zero  lead  time  implies  that  s- =  0. 

Now  item  i,  from  time  to  time,  will  be  faced  with  the  opportunity  of  a 
replenishment  at  the  reduced  setup  cost  of  only  "a-".  Such  an  event  is  caused  by 
another  item  triggering  a  replenishment.  Occasionally  item  i  will  trigger  an  order  when 
its  inventory  hits  the  zero  level.  Under  such  circumstances  we  would  like  to  determine 
the  (S,  c)  pair  which  minimizes  the  expected  costs  per  unit  time. 

Silver  shows  that  the  opportunities  to  replenish  at  the  reduced  cost  occur 
probabilistically  according  to  a  Poisson  process  with  a  rate  ]i  per  year,  where  ]i  is  the 
expected  number  of  orders  triggered  per  year  by  all  other  items  in  the  group. 

Let  NT-  be  the  number  of  orders  triggered  per  year  by  item  i.  Then  \i-  is  given 
by  |X-  =  Y.  NTj,  where  the  summation  is  taken  over  ail  i  #  j. 

For  example,  if  there  are  four  items  in  the  group  and  we  are  considering  item 
2,  we  would  have  ji->  =NT|  +  NT-j  +  NT4. 

It  is  very  difficult  to  determine  the  value  of  p:-  algebraically.  We  will  consider 
some  heuristic  methods  to  approximate  this. 

Silver  [Ref.  3:  eqn.  13]  determines  that  the  relevant  cost  equation  (the 
subscripts  i  have  been  suppressed)  is  given  by: 

Eq=  (S  -  c  +  p  (1  -  pc  )  /  (1  -  p)}-1  {(S  -  c)(S  +  c  +  1)  vr  /  2  (2.20) 

+  p  [c  -  p  (1  -  pc  )  /  (1  -  p)]  vr  /  (1  -  p)  +  X  pc  A  +  X  a} 

where 

p  -  X  /  (X  +  |i) 

The  algorithm  suggested  by  Silver  in  [Ref.  3]  is  a  simultaneous  optimization 
procedure  to  minimize  eqn.  (2.17)  with  respect  to  S-,  c-,  s-.     The  procedure  is  iterative. 

Solution  procedures  to  determine  the  control  parameters  of  the  (S,  c,  s)  policy 
are  explained  in  the  next  chapter. 
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III.  DEVELOPMENT  OF  THE  MODEL 

In  this  chapter  we  consider  the  general  joint  replenishment  problem  with  positive 
leadtimes.     We   obtain  values   of  the   parameters   s-,   c-,   and   S-   using  an  iterative 

procedure  recommended  in  [Ref.  3].  Silver  [Ref.  10]  derives  the  following  expression  for 
the  total  expected  costs  per  year  for  item  i  (we  suppress  the  item  subscript  for 
simplification  of  notation): 

Eq  -  (S  -  c  +  p  (1  -  pc)  /  (1  -  p)}-1  {(S  -  c)(S  +  c  +  1)  vr  /  2  (3.1) 

+  p  [c  -  p  (1  -  pc)/(l  -  p)Jvr/(l  -  p)  +  X  pc  A  +  I  a} 

.  where 

P   =   X;(X  +   ]l) 

p   is   the   probability   that   a   particular  occurance   is   either  a  demand   or   an 
opportunity  to  replenish  at  reduced  cost. 

A.       SOLUTION  PROCEDURE 

We  recommend  as  a  solution  procedure  a  somewhat  modified  version  of  the 
procedure  suggested  by  Silver. 
1.  Algorithm 

Step  1  :  Initialization 

Use  the  values  computed  by  the  independent  control  to  initialize.     Let   k  =  1, 
where  k  =  index  of  iteration. 

cmm  =  °- 1 

cmax  =  Q 

EC(k)  =  EQdetermined  from  the  independent  control) 

N(k)  =  X  /  Q* 

H--(k)  =  ]T  Ny  where  all  i  ¥=  j,  and  i,  j  =  1 n 

p(k)  =  X.  /  (X  +  li) 

c(k)  =  0 

S(k)  =  q" 
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Step  2  •  Let  k  =  k  +  1 

c(k)  =  c(k-l)  +  (cmax  -  crnin)  /  2,  if  EC(k)  <  EC(k-l) 
c(k)  =  c(k-l)  -  (cmax  -  Cmin)  /  2,  if  EC(k)  >  EC(k-l) 
S(k)  =  c(k)  -  p  (1  -  pc),(l  -  p)  +  V  {2  X  (a  +  A  Xc)/vr 

+  2  c  Xc+  l/(l  "  P)  -  P  (1  -P°)0  +  PC+  lW  -  P)2} 
EC(k)  =  {p  (1  -  pc)/(l  -  p)}'1^  [c  -  p  (1  -  pc)/(l  -  p)]vr 
/(  1  -  p)  +  X  pc  A  +  X  a} 

Step  3  :  Test  for  convergence  of  c  and  S. 

Round  the  c  and  S  values  to  the  nearest  integer. 

If  c(k)  =   c(k-l)  and  S(k)  =   S(k-l),  go  to  Step  4,  else  go  to  Step  2.    And, 
repeat  until  S  and  c  converge. 

Step  4  :  Repeat  Steps  1,  2  and  3  for  i  =  1,2,  ...,  n. 

Step  5  :  Let  Nj  =  /.  p  C/[S  -  c  +  p  (1  -  p  c)/(l  -  p)}. 

p  values  are  from  the  previous  iteration  and  S  and  c  values  are  from  step  4. 

Step  6  :  Repeat  steps  1-5  until  |c(k)  -  c(k-l)|<5  and  |S(k)  -  S(k-l)|<5. 

Step  7  :  Now  determine  the  must-buy  points,  s-v  for  i  =    1,  ...,  n  from  the  service 
constraints. 

(a)  For  Service  Measure  1:  probability  (Pi)  of  no  shortage  per  cycle,  find  the 
smallest  value  s  such  that 

(l/p)Pr(X<s  +  dXU-  Pr(X  =  s+l|X  L) 
-ps+1XPr(X  =  x0|XL)(l/p)xo  >  ?vpc 

(b)  For  Service  Measure  2:  fraction  (P->)  of  demand  to  be  satisfied  directly 
from  shelf,  find  the  smallest  value  s  such  that 

pc{U-sXLPr(X<s-l|U)  +  s  Pr(X  <>  s\\  L)} 
+  (1  -  P)c+S  I  P  "wo  (U-w04L  Pr(wQ  -  1|X  L) 
+  w0Pr(X<wQ)} 
<  [S-c+p(l-pc)/(l-p)}(l-P2) 

Step  8  :  Recompute  the  values  S-  and  c-. 

Si  =  S  +  si 

Ci  =  c  +  Si 

Stop 
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2.  Flowchart 


Initial is at ion 

k=l  c-iB-o 

C(k)=0 


k  =  k  *  1 


Sat  trial  of  C 


C(k)=C(k-l)+ 


V 


Computa  S(k)  and  EC(k) 
(Uaa  aqn.a  at  atap  2) 


**C(k) 


CainaC(k) 


C(k) 


I 


T 


Find    pair    of     S,C,W"hen     laad    tima=0 
Cst   ■  C1  )     .  Win .  EC  (  S  j  ,   C  j  ) 


Find    Coordinated    pair    of    5,C 
(Co.puta      NlJit   (  pi    ) 


No 


(      Stop     J 


Figure  3.1     Determination  of  coordinated  pair  (S,c). 
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IV.  SIMULATION  STUDY 

In  order  to  determine  how  costs  of  ordering,  holding,  shortage,  etc.  are  affected 
by  joint  replenishment,  a  simulation  model  was  developed.  The  purpose  of  the 
simulation  study  was  to  test  the  cost  savings  by  comparing  the  costs  for  the  two 
inventory  policies,  independent  and  joint  replenishment.  The  simulation  was  written  in 
FORTRAN77. 

A.  ASSUMPTIONS  FOR  THE  SIMULATION  STUDY 

1.  Poisson  demands  (exponential  interarrival  times). 

2.  Items  demanded  one  at  a  time. 

3.  Constant  lead  times  (independent  of  the  size  of  the  order  and  the  number  of 
items  in  a  replenishment). 

4.  For  the  joint  replenishment  case,  place  an  order  whenever  anv  item  hits  its 
must-buv  point  s-  and  order  up  to  S-.  Also,  include  item  j  #  i  if  the  inventorv 
position'for  item]  is  below  the  can- buy  point  c-.   Order  up  to  S-. 

5.  For  the  independent  control,  order  up  to  level  S:  whenever  the  stock  on  hand 
for  item  i  reaches  s-. 

6.  For  a  given  item,  there  is  never  more  than  a  single  order  outstanding  in  the 
joint  replenishment  case. 

B.  MODEL  INPUTS 

1.  Number  of  items. 

2.  Demand  rate  "k-  for  each  item. 

3.  Holding  cost  rate  v-  for  each  item. 

4.  Lead  time  (in  months)  L. 

5.  Type  of  service  level  constraint  (Typel,  or  Type2) 

6.  Service  level  required. 

7.  Unit  price  UPj. 

8.  Group  ordering  cost  A. 

9.  Individual  ordering  cost  a-. 

10.  Stockout  cost  rate  per  unit  stockout  (percentage  of  unit  price). 

11.  Time  weighted  stockout  cost  (percentage  of  unit  price). 

12.  S-,  c-,  s-  for  each  item. 

13.  VIodel  options  (Independent  control  or  Joint  Replenishment,  etc.) 

14.  Limits  of  simulation  period,  NSIM. 
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C.  MODEL  OUTPUTS 

1.  Achieved  service  level  for  each  item. 

2.  Number  of  orders  (independent  and  joint). 

3.  Total  ordering  cost  for  each  item. 

4.  Total  holding  cost  for  each  item. 

5.  Total  unit  years  of  backorders. 

6.  Total  number  of  backorders. 

7.  Total  unit  years  of  stock  held. 

8.  Total  stockout  cost  (fixed  and  time  weighted)  for  each  item. 

9.  Total  annual  costs  for  each  item. 

10.  Total  annual  variable  costs  for  all  items. 

11.  Standard  deviation  of  total  annual  variable  costs. 

12.  A  complete  audit  trail  of  all  orders  placed  (optional). 

D.  STRUCTURE  OF  THE  PROGRAM 

1.  Algorithm 

Step  1  :  Read  input  data,  process  options  and  initialize  variables. 

Step  2  :  Determine  the  earliest  event  time  and  type  of  event. 

Let  TDA-  be  the  time  of  the  next  demand  for  item  i  and  TOA-  be  the  next 
time  of  the  for  item  of  an  order  for  item  i. 

Set  the  master  clock  time,  CT,  to  the  minimum  of  {TDA^  TOA-}. 

If  the  event  type  is  a  demand  for  item  i  then  generate  another  demand  time 
for  that  item. 

If  the  event  type  is  an  order  arrival  then  increament  the  on  hand  inventory  for 
all  items  included  in  the  order. 

Step  3  :  Determine  if  an  order  is  required. 

If  the  onhand  level  is  less  than  or  equal  to  the  reorder  point  for  that  item  and 
there  are  no  outstanding  orders  for  that  item,  place  an  order.  If  the  ordering  policy  is 
independent,  go  to  step  5. 

Step  4  :  Test  for  joint  replenishments. 

If  the  onhand  level  is  less  than  or  equal  to  the  can  order  point  for  another 
item  j  in  the  replenishment  set,  place  joint  replenishment  for  the  items  to  the  level  S-. 
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Step  5  :  Determine  if  the  simulation  is  complete. 
If  CT  ^  NSIM  stop,  else  go  to  step  2. 

The  program  collects  detailed  statistics  at  each  event  time  for  unit  years  of 
stock,  held,  number  of  stockouts,  unit  years  of  stockouts,  number  of  orders,  number  of 
units  ordered,  and  number  of  demands  for  each  item. 

At  the  conclusion  of  the  simulation  additional  summary  statistics  are 
computed  for  achieved  service  levels  and  costs. 
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2.  Flowchart 
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Figure  4.1     Simulation  Flowchart. 
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V.  EVALUATION  OF  THE  MODEL 

A.       PERFORMANCE  ANALYSIS 
1.  Basic  Examples  Tested 

Table  4  shows  the  data  used  for  the  evaluation  of  the  joint  replenishment 
model.  Several  combinations  (96  cases)  of  service  level,  service  measure,  and  lead  times 
were  considered.     The  data  include  the  test  set  considered  in  Silver's  paper  [Ref.  3]. 

TABLE  4 
ITEM  CHARACTERISTICS 


group 

A 

a 

number 

ITEM  DATA 

number 

(S) 

(S) 

of  items 

l 

/.j  (piece/yr.)       Vj 

(S/  piece) 

1 

50 

10 

4 

1 

2 

i 

-t 

290.000 
41.000 
77.000 

122.000 

6.90 
1.20 
3.90 
2.30 

2 

125 

so 

10 

1 
2 

3 

4 
5 
6 
7 
8 
9 
10 

40.630 
4.090 

34.680 
4.240 
4.240 
4.090 
4.240 

28.780 
4.090 
4.090 

0.45 

36.41 

42.17 

4.46 

4.46 

36.41 

4.46 

3.92 

36.41 

36.41 

3 

50 

5 

4 

1 

j 
4 

290.000 
41.000 
77.000 

122.000 

6.90 
1.20 
3.90 
2.30 

4 

50 

5 

8 

1 

2 

4 
5 
6 

7 
8 

290.000 
41.000 
77.000 

122.000 
50.000 

154.000 
87.000 
25.000 

6.90 
1.20 
3.90 
2.30 
1.20 
3.90 
2.30 
1.20 

2.  Model  Results 

a.   Convergence  of  Parameter  Values  and  Costs 

Table  5  shows  the  results  of  each  iteration  from  Silver's  model  for  joint 
replenishment.     Observe    that   live   iterations   were   required   for   convergence.      This 
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example  considers  service  measure  type  2  with  a  specified  level  0.99.     The  lead  time  is 
one  month. 

TABLE  5 
PARAMETERS  AND  COSTS 


Item 

Iter. 

P 

C 

S 

EC 

1 

1 

0.9942 

79 

146 

201.8498 

2 

0.9972 

84 

152 

210.4767 

3 

0.9973 

84 

152 

210.5555 

4 

0.9974 

87 

152 

210.8846 

5 

0.9973 

87 

152 

210.8239 

2 

1 

0.9272 

31 

90 

21.7035 

2 

0.9499 

36 

97 

23.4369 

i 
j 

0.9464 

36 

96 

23.1011 

4 

0.9469 

36 

96 

23.1429 

5 

0.9467 

36 

96 

23.1271 

3 

1 

0.9650 

27 

85 

66.6963 

• 

2 

0.9768 

48 

89 

70.0377 

-* 
j 

0.9741 

43 

88 

68.S164 

4 

0.9744 

43 

88 

68.9584 

5 

0.9744 

43 

88 

68.9309 

4 

1 

0.9774 

45 

138 

63.6090 

2 

0.9850 

71 

145 

66.8895 

-> 

j 

0.9835 

70 

143 

65.8678 

4 

0.9835 

70 

143 

65.9083 

5 

0.9835 

70 

143 

65.8689 

b.    Cost  Savings'.  Joint  Replenishment  VS.  Independent  Control 

Table  6  shows  the  cost  savings  achieved  by  the  joint  replenishment  policy 
for  both  the  modified  Silver's  method  and  our  heuristic  method. 
The  comparisions  are  provided  for  two  cases: 

1.  Service  measure  1,  Pj  =  0.95 

2.  Service  measure  2,  P->  =  0.99 

Observe  that  both  the  heuristic  method  and  Silver's  method  achieve  cost 
savings  greater  than   13  %  by  using  joint  replenishment,  and  the  heuristic  method 
results  in  greater  cost  savings  than  does  Silver's  method. 
3.  Simulation  Results 

Figure  5.1  shows  the  time  history  of  stock  on  hand  for  items  1  to  4  generated 
by  one  simulation.  The  upper  graph  represents  individual  control  and  the  lower  graph 
represents  joint  replenishment. 


30 


TABLE  6 
COMPARISION  OF  TWO  METHODS 


1.  modified  Si  J 

.ver '  s 

method  : 

pi  =  o- 

95,  L  =  1  month 

item 

P 

s 

C 

s 

EC( JP) 

EC( IP) 

Save(%) 

1 

0. 9973 

32 

119 

184 

221. 5339 

232. 024 

4.  48 

2 

0. 9467 

4 

40 

100 

23. 2571 

35. 343 

34.  17 

3 

0. 9744 

9 

52 

97 

70. 9459 

88. 860 

20.  16 

4 

0. 9835 

13 

83 

156 

67. 1722 

84.  977 

20.  95 

total  cost 

383. 0190 

441. 203 

13.  19 

2.  heuristic  me 

ithod  : 

pl  =  o- 

95,  L  = 

1  month 

item 

P 

s 

c 

S 

EC( JP) 

EC( I?) 

Save(%) 

1 

0. 9942 

31 

110 

■  177 

211. 2798 

232. 024 

3.  94 

2 

0. 9272 

3 

34 

93 

21.  6035 

35. 3  43 

3.3.  37 

3 

0. 9650 

9 

36 

94 

58. 7113 

38. 360 

22.  57 

4 

0. 9774 

14 

59 

152 

65. 3723 

84. 977 

23.  07 

tOtc 

il  cost 

351. 3467 

441. 203 

16.  33 

3.  modified  Si] 

.ver' s 

method  : 

P2  =  o. 

99,  L  =  1  month 

item 

P 

s 

C 

s 

EC(  J? ) 

EC( IP) 

Save(%) 

1 

0. 9973 

25 

112 

177 

211. 9739 

220. 984 

4.  08 

2 

0. 9467 

-1 

35 

95 

22. 0671 

34. 383 

35.  82 

3 

0. 9744 

5 

48 

93 

67. 8259 

85. 740 

20.  39 

4 

0. 9835 

7 

77 

150 

64. 4122 

32. 217 

21.  56 

totz 

il  cost 

365. 2788 

423. 323 

13.  48 

4.  heuristic  me 

ithod  : 

P2  =  o. 

99,  L  = 

1  month 

item 

P 

s 

c 

S 

EC( JP) 

EC( I?) 

Save(%) 

1 

0. 9942 

25 

104 

171 

202. 9998 

220. 984 

8.  14 

2 

0. 9272 

-1 

30 

89 

20. 6435 

34.  383 

39.  96 

3 

0. 9650 

5 

32 

90 

65. 5912 

85. 740 

23.  50 

4 

0. 9774 

8 

53 

146 

62. 6123 

82. 217 

23.  35 

tots 

il  COSt 

351. 8467 

423. 323 

16.  88 
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The  simulation  results  for  several  combinations  of  leadtime.  service  measure 
and  service  level  are  represented  in  Tables  7,8,9. 

When  lead  time  is  less  than  1  month,  the  heuristic  method  yields  more  savings 
in  total  cost  and  also  satisfies  the  required  service  level.  But,  as  the  lead  time  becomes 
larger,  it  sometimes  fails  to  meet  the  required  service  level,  although  it  does  continue  to 
generate  greater  cost  savings  than  Silver's  method  with  the  cost  parameters  examined 
in  this  thesis.  Both  methods  generate  substantial  savings  over  independent  control. 
In  short  lead  times,  both  methods  tend  to  overprotect  against  stockouts.  They 
generally  produce  service  levels  in  excess  of  what  is  required. 

One  explanation  for  the  better  performance  of  the  joint  replenishment  models 
for  the  short  lead  time  examples  is  the  fact  that  a  key  step  in  the  joint  replenishment 
algorithm  for  determining  the  S,  c,  s  parameters  assumes  that  lead  time  is  zero. 
Adjustments  are  then  made  to  those  values  to  account  for  the  nonzero  lead  times  but, 
apparently  those  adjustments  are  not  adequate.  Additional  work  is  required  to 
examine  this  issue  and  to  explore  other  methods  of  adjustment. 

We  should  keep  in  mind  the  fact  that  the  model  satisfies  all  requirements  in 
the  experiment  does  not  mean  that  it  is  a  correct  model  —  merely  that  it  is  a  plausible 
one  which  has  not  been  found  inadequate  by  the  data  or  experiment. 
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2  YEARS  RUN  :  ITEMS  IN  GROUP   1 
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Figure  5.1     Stock  on  Hand  —  Independent  Control  and  Joint  Replenishment. 
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TABLE  7 
COST  SAVINGS  BY  JOINT  REPLENISHMENT:  WHEN  L  =  1  MONTH 


1.     PARAMETERS    FROM   MODIFIED    SILVER'S    METHOD 


Group 
number 

req. 

SVC 

(%) 

1 

Pl=95 
Pl=99 

P2=95 
P2=99 

2 

Pl=95 
Pl=99 

P2=95 
P2=99 

3 

Pl=95 
Pl=99 

P2=95 
?2=99 

4 

Pl=9  5 
Pl=99 

P2=95 
P2=99 

AVG.    Total   Variable   C0sts($) 


INDEPENDENT 


MODEL 


441. 20 
448. 06 


403. 20 
423. 32 


1111. 43 
1153. 33 


1032. 91 
1088. 91 


423. 30 
430. 16 


386. 08 
406. 30 


676. 10 
686. 94 


620. 56 
650. 70 


SIM 

443. 

445. 

04 

77 

414. 
431. 

71 
00 

1073. 
1107. 

85 
24 

1029. 
1066. 

70 
90 

422. 
425. 

05 

74 

396. 
409. 

72 
34 

677. 
683. 

55 
89 

627. 
650. 

49 
30 

JOINT 

MODEL 

SIM 

383. 
392. 

02 

74 

361. 
368. 

48 
70 

343. 
366. 

31 
23 

341. 
348. 

12 

38 

998. 
1037. 

32 
70 

953. 
956. 

78 
40 

878. 
971. 

37 
82 

888. 
926. 

82 
96 

356. 
365. 

89 
S3 

345. 
354. 

93 

15 

315. 
340. 

71 

01 

315. 
323. 

35 

05 

513. 
528. 

14 
48 

490. 
508. 

92 
80 

453. 
490. 

34 
60 

433. 
459. 

05 
33 

Cost  Savings(%) 
by  Joint 
Replenishment 

MODEL 

SIM 

13. 
12. 

19 

35 

13.  41 
17.  29 

14. 
13. 

86 

48 

17.  75 
19.  17 

10. 
10. 

18 
03 

11.  22 
13.  62 

14. 
10. 

96 
75 

13.  63 
13.  12 

15. 
•  14. 

69 
95 

18.  04 
16.  32 

18. 
16. 

23 
42 

20.  51 

18.  74 

24. 
23. 

10 
07 

27.  54 
25.  60 

26. 
24. 

14 
60 

30.  99 
29.  42 

AVG.  SVL 

achieved 

bv 
SIM. 

99. 
99. 

49 
74 

97. 
99. 

91 
71 

99. 
100. 

27 
00 

97. 
99. 

70 
64 

99. 
100. 

76 
00 

98. 
99. 

44 

76 

99. 

99. 

21 
30 

99. 
99. 

01 
33 

2.     PARAMETERS    FROM   HEURISTIC    METHOD 


Pl= 
Pl= 

=95 

=99 

P2  = 
P2  = 

=95 
=99 

PI: 

PI: 

=95 
=99 

P2  = 
P2  = 

=95 
=99 

Pl: 
Pl= 

=95 
=99 

P2  = 

P2: 

=95 
=99 

Pl  = 
PI: 

=95 
=99 

P2  = 

P2  = 

=95 
=99 

441. 
443. 

20 
06 

443. 
445. 

04 
77 

403. 
423. 

20 
32 

414. 
431. 

71 
02 

1111. 
1153. 

43 
33 

1073. 
1107. 

35 

24 

1032. 
1088. 

91 
91 

1029. 
1066. 

70 
90 

423. 
430. 

30 

16 

422. 
425. 

05 

74 

386. 

406. 

08 
80 

396. 
409. 

72 

84 

676. 
686. 

10 
94 

677. 
633. 

55 
89 

620. 
650. 

56 
70 

627. 
650. 

49 
80 

366. 
376. 

97 
93 

329. 
351. 

23 
85 

949. 
999. 

18 
68 

849. 
903. 

23 

69 

332. 
343. 

99 
17 

296. 
319. 

95 

11 

447. 
463. 

50 
88 

392. 
426. 

74 
24 

365. 30 
363. 73 


329. 
347. 

36 
15 

960. 
979. 

01 

11 

845. 
893. 

07 
75 

342. 
350. 

74 
13 

314. 
323. 

70 
37 

455. 
471. 

68 
45 

394.  17 
420. 22 


16. 
15. 

S3 
88 

18. 
16. 

35 
88 

14. 
13. 

60 
32 

17. 
17. 

78 
01 

21. 
20. 

33 
22 

23. 
21. 

09 
56 

3"3 
32! 

81 
47 

36.  71 
34.  50 


17.  55 
17.  28 


20. 
19. 

53 
46 

10. 
11. 

60 

57 

17. 
16. 

93 
23 

18. 
17. 

79 
76 

20. 

19. 

67 
76 

32. 
31. 

75 
06 

37.  18 
35.  43 


96.  51 
99.  04 


97. 
99. 

34 
46 

97. 
99. 

43 
44 

97. 
99. 

54 
08 

93. 
99. 

60 
77 

98. 
99. 

01 

64 

97. 
99. 

19 
51 

97.  54 
99.  29 
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TABLE  8 
COST  SAVINGS  BY  JOINT  REPLENISHMENT:  WHEN  L  =  3  MONTH 


1.  PARAMETERS  FROM  MODIFIED  SILVER'S  METHOD 


Group 
number 

req. 
SVL 
(%) 

1 

Pl=95 
Pl=99 

P2=95 
P2=99 

2 

Pl=95 
Pl=99 

P2=95 
P2=99 

3 

Pl=95 
Pl=99 

P2=95 
P2=99 

4 

Pl=95 
Pl=99 

P2=95 
P2=99 

AVG.  Total  Variable  COsts($) 


INDEPENDENT 


MODEL 


453. 
466. 

70 
64 

404. 

431. 

20 
48 

1173. 
1224. 

94 
26 

1047. 
1109. 

78 

34 

435. 
•  448. 

80 

.74 

386. 

413. 

76 
82 

694. 
713. 

66 
84 

619. 70 
660. 82 


SIM 

451. 
464. 

55 
38 

407. 
433. 

47 
39 

1135. 
1177. 

01 
53 

1054. 
1097. 

40 
29 

439. 
446. 

10 
89 

390. 
410. 

60 
99 

690. 
711. 

92 
89 

623. 99 
660. 38 


JOINT 


MODEL 


394.  04 
408. 46 


340. 27 

373. 06 


1019. 79 
1078. 73 


893. 13 
992. 50 


365. 75 
382. 33 


314.  71 
347. 03 


525. 31 
549. 89 


446. 77 
498. 85 


SIM 

368. 
383. 

34 
34 

331. 
350. 

60 
99 

938. 
994. 

40 

54 

855. 
917. 

98 
02 

358. 
371. 

36 
04 

309. 
337. 

32 

40 

491. 
516. 

81 
17 

418. 
467. 

44 
03 

Cost  Savings(%) 
by  Joint 
Replenishment 

MODEL 

SIM 

13.  15 
12.  47 

18.  43 

17.  45 

15.  82 
13.  54 

18.  62 

19.  01 

13.  13 
11.  89 

17.  32 
15.  54 

14.  76 
10.  53 

18.  32 
16.  42 

16.  07 
.  14.  80 

13.  39 
16.  97 

13.  63 
16.  14 

20.  31 

17.  91 

24.  38 
22.  97 

23.  32 
27.  49 

27  gi 

24!  51 

32.  94 
29.  28 

AVG. SVL 

achieved 

bv 

SIM. 


98. 
99. 

74 

96. 
99. 

91 
59 

98. 
99. 

33 
33 

98. 
99. 

24 
73 

96. 

99. 

73 
51 

97. 
99. 

38 

57 

99. 
99. 

32 
80 

98.  70 

99.  84 


PARAMETERS  FROM  HEURISTIC  METHOD 


PI: 
Pl  = 

=95 
=99 

P2  = 
P2  = 

=95 
=99 

Pl  = 
PI: 

=95 

=99 

p2: 

P2: 

=95 
=99 

PI: 
PI: 

=9  5 
=99 

P2  = 
P2  = 

=95 
=99 

PI: 
Pl  = 

=95 
=99 

P2  = 
P2  = 

=95 
=99 

453. 70 
466. 64 

404.  20 
431.  48 

1173. 94 
1224.  36 


1047. 78 
1109. 34 

435. 30 
448. 74 

386. 76 
413. 82 

694. 66 
713. 84 

619. 70 
660. 82 


451. 55 
464. 38 


407. 47 
433. 39 


1135. 01 
1177. 53 


1054. 40 
1097. 29 

439. 10 
446. 39 


390. 60 
410. 99 


690. 92 
711. 89 


623. 99 
660. 38 


377. 
392. 

98 
64 

326. 
357. 

78 
32 

973. 
1041. 

83 

41 

354. 
950. 

13 
04 

340. 
353. 

91 
43 

291. 
324. 

57 

41 

456. 
483. 

20 
56 

381. 
432. 

62 
52 

363. 73 
378. 46 


318. 
346. 


15 
98 


966. 47 
1015. 42 

360. 07 
923. 59 


3  44 
360 


09 
38 


3  01. 49 
327. 00 

447. 29 
474. 04 

400. 09 
425. 76 


16. 
15. 

69 
36 

19. 
17. 

15 
19 

16. 
14. 

52 
94 

18. 
14. 

48 
36 

21. 
20. 

77 

13 

24. 
21. 

61 
61 

34. 
32. 

33 

26 

38.  42 
34.  55 


19.  45 
IS.  50 


21.  92 
19.  94 


14.  35 
13.  77 


18.  43 
15.  83 


21.  64 
19.  25 


n  go 
20!  44 


35.  26 
33.  41 


35.  88 
35.  53 


97. 
99. 

35 

76 

96. 
99. 

75 
44 

*93 
99'. 

78 
17 

97. 
99. 

94 
44 

99. 
99. 

30 
77 

97. 

99. 

70 
74 

95. 
98. 

34 
97 

95.  55 
99.  10 
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TABLE  9 
COST  SAVINGS  BY  JOINT  REPLENISHMENT:  WHEN  L  =  6  MONTH 


1.     PARAMETERS    FROM   MODIFIED    SILVER'S   METHOD 


Group 

number 

req. 
SVL 
(%) 

1 

Pl=95 
Pl=99 

P2=95 
P2=99 

2 

Pl=95 
Pl=99 

P2=95 
P2=99 

3 

Pl=95 
Pl=99 

P2=95 
P2=99 

4 

Pl=95 
Pl=99 

P2=95 
P2=99 

AVG.     Total   Variable   C0sts($) 


INDEPENDENT 


MODEL 

SIM 

465. 
485. 

47 
50 

461. 
482. 

92 
24 

406. 

25 

41 

407. 
440. 

04 
34 

1217. 
1285. 

73 
89 

1174. 
1240. 

03 
05 

1046. 
1161. 

00 
47 

1056. 
1132. 

59 

14 

A47 
467.' 

56 

.60 

442. 
464. 

91 

06 

390. 

423. 

96 
50 

392. 
418. 

25 
99 

712. 
740. 

21 
67 

711. 
737. 

50 
77 

624. 
674. 

71 
81 

629. 
674. 

18 
72 

JOINT 


MODEL 

SIM 

404. 
425. 

32 

84 

374. 
396. 

87 

24 

343. 
381. 

22 
50 

320. 
354. 

37 
81 

1063. 
1131. 

46 
83 

983. 
1045. 

25 

79 

897. 
1014. 

43 
36 

851. 
933. 

39 

91 

375. 
397. 

58 
10 

359. 
■  380. 

91 
81 

317. 
355. 

44 
48 

310. 
3  42. 

39 
46 

536. 
571. 

58 

58 

505. 
529. 

25 
95 

4.4.7 
5"d7! 

04 
72 

411. 
465. 

24 
69 

Cost  Savings(%) 
by  Joint 

Rep leni  shment 

MODEL 

SIM 

13. 
12. 

14 
29 

18.  35 
17.  83 

15. 

13. 

51 
57 

21.  29 
19.  42 

12. 
11. 

26 
98 

16.  25 
15.  67 

14. 
12. 

20 
67 

19.  42 
17.  51 

15. 
15. 

08 
08 

18.  74 
17.  94 

18. 
16. 

31 
06 

20.  87 
18.  27 

24. 
22. 

66 
83 

28.  99 
28.  17 

23. 
24. 

44 
76 

34.  64 
30.  98 

AVG.  SVL 

achieved 

bv 
SIM. 

97. 

99. 

53 
24 

96. 
99. 

23 
34 

97. 
99. 

90 
33 

98. 
99. 

07 
77 

94. 
99. 

30 

14 

97. 
99. 

19 
29 

96. 

99. 

47 
04 

97. 
99. 

19 

44 

2.  PARAMETERS  FROM  HEURISTIC  METHOD 


Pl= 
Pl= 

=95 
=99 

P2  = 
P2= 

=95 
=99 

Pl= 

Pl= 

=9  5 
=99 

P2  = 

P2  = 

=95 
=99 

Pl  = 

Pl  = 

=95 
=99 

P2  = 

P2: 

=95 
=99 

Pl  = 
Pl= 

=95 
=99 

P2  = 
P2  = 

=95 
=99 

465. 47 
485. 50 


406.25 

441. 41 

1217. 73 
1285. 89 


1046. 00 
1161. 47 

447. 56 
467. 60 

390. 96 
423. 50 

712. 21 
740. 67 

624. 71 
674.  81 


461. 92 
482. 24 

387. 43 
409. 79 

407. 04 
440. 34 

328. 41 
366. 23 

1174. 03 
1240. 05 

1030. 96 
1090. 96 

1056. 59 
1132. 14 

848. 43 
962. 32 

442. 91 
464. 06 

351. 19 
372. 95 

392. 25 

418. 99 

292. 41 

331.  47 

711. 50 
737. 77 

466. 69 
503. 17 

629. 18 
674.  72 

381.  93 
442. 41 

379. 63 
401. 59 


328. 97 
363. 02 


991. 45 
1049. 95 

878. 96 
940. 00 


341. 63 
362. 02 


298. 23 
330. 47 


532. 09 
559. 74 


404. 55 
452. 60 


16.  77 
15.  60 


19. 
17. 

16 
03 

15. 
15. 

34 

16 

IS. 
17. 

39 
15 

21. 
20. 

53 
24 

25. 

21. 

21 
73 

34. 

32. 

47 
07 

38.  36 
34.  44 


17.  32 
16.  72 


19. 
17. 

18 
56 

15. 
15. 

55 
33 

16. 
16. 

81 
97 

22. 

21. 

37 
99 

23. 

21. 

97 
13 

25. 

24. 

21 

13 

35.  70 
32.  92 


*92. 64 
98.  74 


95. 
99. 

33 

08 

93. 
100. 

92 
00 

96. 
99. 

53 
34 

*91. 

*96. 

90 

14 

95. 
93. 

57 
48 

*86. 
*94. 

62 
34 

*92. 36 

*97. 70 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

For  each  of  several  test  combinations,  we  determined  the  independent  and 
coordinated  control  parameter  values  and  the  associated  expected  costs.  We  simulated 
several  years  of  operations  with  the  derived  joint  replenishment  ordering  parameters  to 
evaluate  the  effectiveness  of  the  models  as  compared  to  independent  control  models. 
We  also  use  the  simulation  to  compare  a  simple  heuristic  joint  replenishment  model. 
With  the  simulation  we  were  able  to  accumulate  costs  (ordering,  holding  and  stockout) 
and  effectiveness  as  measured  by  two  service  levels. 

In  summary,  the  following  statements  can  be  made: 

1.  Coordinated  control  results  in  total  costs  significantly  lower  than  under 
independent  control.   The  average  cost  savings  are  about  T8.7  %  . 

2.  Coordinated  control  provides  better  service  than  independent  control.  In  a  few 
cases,  the  actual  service  for  an  individual  item  is  sliehtlv  lower  than  required: 
when  this  occurs,  the  item  involved  is  the  one  that  triggers  the  majority  of  the 
replenishments  in  the  group. 

3.  Cost  savings  of  coordinated  control  increases  as  the  ratio  a/ A  decreas    es. 

4.  Cost  savings  improve  as  the  number  of  items  in  replenishment  set  increases:  the 
more  items  there  are  in  the  group,  the  more  beneficial  coordinated 
replenishment  becomes. 

5.  For  a  given  service  measure  and  lead  time,  the  cost  savings  increase  as  the 
desired  service  level  decreases. 

B.  RECOMMENDATIONS 

The  procedure  to  determine  the  parameters  for  coordinated  control  requires  only 
the  data  for  independent  control.  We  observe  that  the  joint  replenishment  models 
perform  best  for  shorter  lead  times.  This  is  probably  the  result  of  the  fact  that  the 
joint  replenishment  algorithms  initially  assume  zero  lead  times  and  then  adjust  the 
derived  parameters  heuristically  to  account  for  nonzero  lead  times.  This  adjustment 
procedure  may  be  inadequate.  Additional  work,  is  needed  to  consider  other  adjustment 
methods. 
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APPENDIX  A 
"PARA":PROGRAM  TO  SELECT  PARAMETERS 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


** 


VARIABLE  DEFINITION   ** 


11 


DM 

MR 

S/Sl 

C/Cl 

IPC 

IPS 

ECI 

EC/EC1 

TLM 

TLT 

TL 

DLT 

NIT 

TN 

R  / 

F 

UP 

OCI 

OCG 

SVL 

SCU 

PSO 

STON 


RHO 


ANNUAL  DEMAND 

REORDER  POINT 

ORDER  -  UP  POINT 

CAN  ORDER  POINT 

OPTIMAL  PAIR  OF  C 

OPTIMAL  PAIR  OF  S 

ECNOMIC  ORDER  QUANTITY 

AVERAGE  ANNUAL  TOTAL  VARIABLE  COST (INDEPENDENT  ODER) 

AVERAGE  ANNUAL  TOTAL  VARIABLE  COST (JOINT  ORDER) 

LEAD  TIME (IN  MONTH) 

LEAD  TIME (IN  YEARS) 

AVERAGE  LEAD  TIME (IN  YEAR) 

LEAD  TIME  DEMAND 

NO.  OF  TOTAL  ITEM  IN  THE  GROUP 

NO.  OF  TRIGGERING  AN  ORDER  BY  ITEM  I  IN  THE  GROUP 

OPPORTUNITY  TO  REPLENISH  AT  REDUCED  COST 

HOLDING  COST  RATE  TO  PROCUREMENT  COST 

UNIT  PRICE 

INDIVIDUAL  ORDERING  COST 

GROUP   ORDERING  COST 

SERVICE  LEVEL 

IMPUTED  UNIT  STOCKOUT  COST  GIVEN  SVL. 

MAX.  ALLOWABLE  ROBABILITY  OF  STOCK  OUT 

EXPECTED  NUMBER  OF  STOCKOUT  DURING  A  ORDER  CYCLE 


DIMENSION  ITM(  100 ),  DM(  100 ),  UP  ( 100 ),  H(  100  )  ,OCI(  100  ),ECI(  100),  0.(100) 
&,C(100),TN(100) .S(100),IS(100),R(100),U(100),EC(100),MR(100), 
&IQ ( 100 ), DLT ( 100 ),SCU( 100) ,10(100), RHO (50 ,100), ECI (50, 100), 
&C1(50,100),51(50,100),IC1(50/100),IS1(50,100),IPC(100),IPS(100), 
&MR1(100),TL(100) 

*  DATA  INITIALIZATION 
STN=0 . 

TVC=0 . 

*  READ  PROCESS  OPTIONS 

CALL  OPTION ( ID , LT , KEY , LTT , SVL , MSRCH) 
READ(ID,110)NIT,OCG,F 
TLT=LT/12. 
DO  11  1=1, NIT 

READ(ID,111)ITM(I),DM(I),UP(I),0CI(I) 
t      H(I)=UP(I)*F 
IF(LTT.EQ.1)TL(I)=TLT 
IF(LTT.EQ.2)THEN 
CALL  LEAD(TLT,TLY) 
TL(I)=TLY 
END  IF 
CONTINUE 

*  PRINT  SELECTED  OPTIONS 

WRITE (7,100) KEY , SVL , LTT , ID , MSRCH 

*  CALCULATE  E.O.Q 

DO  1  I  =1,NIT 
DLT(I)=DM(I)*TL(I) 
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OC=OCG+OCI(I) 

o(i)»s6rt72*6m(i)*oc/h(i)) 
iq(i)=int(q(i)+.5) 
c 

C  *   COMPUTE  REORDER  POINT :IND.  ORDER  POLICY 

C 

C  KEY  =  1  :  PF.CE.OF  NO  SHORTAGE  PER  REPLENISHMENT  CYCLE 

C  KEY  =  2  :  FRACTION  OF  DEMAND  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF 


C 


CALL  POISON(KEY,DLT(I),Q(I),SVL/MX/M/MAX,CCDFl) 

MR1(I)=M 
IF(KEY.EQ.1)THEN 

SCU(I)=H(I)*Q(I)/((1.-  SVL)*DM(I)) 
END  IF 
IF(KEY.EQ.2)THEN 

SCU(I)=H(I)*Q(I)/(CCDF1*DM(I)) 
END  IF 

ECI(I)=H(I)*((Q(I)+l.)/2.+(MRl(I)-DLT(I)))+OC*(DM(I)/Q(I)) 
WRITE(7,501)ITM(I),DLT(I),TL(I),Q(I)/IQ(I)/MR1(I),ECI(I),SCU(I) 

TVC=TVC+ECI(I) 


TN(I)=DM(I)/Q(I) 
STN=STN+TN(I) 

1  CONTINUE 
C 

WRITE(7,502)TVC 
C 

C    ■  *   DETERMINE  THE  PARAMETERS  FOR  JOINT  ORDER  POLICY 
C  (REORDER  /  CAN  ORDER  /  ORDER-UP  POINT) 

C 

DO  2  1=1, NIT 

2  U(I)=STN-TN(I) 
DO  3  1=1, NIT 

3  R(I)=DM(I)/(DM(I)+U(I)) 
C 

C     *  SEARCH  LOOP  FOR  DETERMINIG  A  PAIR  (S,C)  :  ITEM  BY  ITEM 
C 

11=1 
SEC=0. 
SEC1=0. 
41    DO  4   1=1, NIT 

RHO(I,II)=R(I) 
K=l 

CMIN=0 . 
CMAX=Q(I) 
C(K)=0. 
S(K)=Q(I) 
EC(K)=ECI(I) 
TEST=10.**10 
ICASE=1 
LOOP=0 
5   K=K+1 

IF(ICASE.EQ.1)THEN 

C(K)=C(K-l)+(CMAX-CMIN)/2 
END  IF 
IF(ICASE.EQ.2)THEN 

C(K)=C(K-l)-(CMAX-CMIN)/2 
END  IF 

S(K)  =  C(K)-R(I)*(l-R(I)**C(K))/(l-R(I))+SQRT(2*DM(I)*(OCI(I) 
&  +OCG*R(I)**C(K))/H(I)+2*CrK)'tR(I)**(C(K)  +  l)/(l-R(I)' 

&  -R(I)*(1-R(I)**C(K))*(1+R(I  **(C(K)+1))/(1-R(I))**2' 


EC  K  =  1./  S  K  -C (K  +R  I  *  l.-R (I  **C (K)  /  l.-R  I  ) 
&*((S(K)-C(K))*(S(K)+C(K)+l)*H(I)/2+R(I)*(C(K)-R(I)*(l.-R(I)**C(K)) 
&/(1.-R(I)))*H(I)/(1.-R(I))+DM(I)*R(I)*Mk)*OCG+DM(I)*OCI(I)) 

IS(K)=INT(S(K)+.5) 

IC  K)=INT(C(K)+.5) 

IF(C(K) i.EQ.C(K-l) .AND.S(K) .EQ.S(K-l) )GO  TO  6 

IF(EC(K) .GT. (0.9995*EC(K-1)) .AND.EC(K) .LE.EC(K-l))  GO  TO  6 

IF(EC(K).LT.EC(K-1))  THEN 

IF(C(K).GT.C(K-1))THEN 
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CHIN=C(K) 

ICASE=1 
GO  TO  5 
ELSE 

CMAX=C(K) 
ICASE=2 
GO  TO  5 
END  IF 
ELSE 

IF(C(K) .LE.C(K-1))THEN 
CMIN=C(K) 
ICASE=1 
GO  TO  5 
ELSE 

CMAX=C(K) 
ICASE=2 
GO  TO  5 
END  IF 
END  IF 
6    IF(LOOP.EQ.0)THEN 
CMIN=0 
CMAX=C ( 2 ) 
ICASE=2 
LOOP=l 
GO  TO  5 
END  IF 
C 
C    *  FIND  OUT  THE  PAIR  (S,C),  MINIMIZES  EC 


C 


C 


DO  7  J=1,K 

IF(EC(J).LE.TEST)THEN 
TEST=EC(J) 
EC1(I,II)=EC(J) 
C1(I,II)=C(J) 
S1(I,II)=S(J) 
END  IF 
CONTINUE 

ICl(I,II)=INT(Cl(I,II)+.5 
ISl(I,II)=INT(Sl(I,II)+.5 
SEC1=SEC1+EC1(I,II) 
4     CONTINUE 

SHAVE=(TVC-SEC1)/TVC 

WRITE(7, 507)11 
DO  20  1=1, NIT 

WRITE(7/508)I,RHO(I,II)/IC1(I/II)/IS1(I,II)/EC1(I,II) 
20    CONTINUE 

WRITE (7,509) SEC1 , SHAVE 
C 

IF(MSRCH.EQ.2)  GO  TO  42 
C 
C    *   COMPUTATIONS  NEEDED  ONLY  FOR  S.  METHOD  1 


C 


c 


IF(II.GE.2)THEN 
DO  21  1=1, NIT 

21  IF(C1(I,II).EQ.C1(I,II-1).AND.S1(I,II).EQ.S1(I,II-1))  GO  TO  42 
END  IF 

IF(II.GE.2)THEN 
DO  22  1=1, NIT 

22  IF(II.GT.NIT.0R.EC1(I,II).GT.(0.9995*EC1(I,II-1)) 

&  .AND.EC1(I,II).LE.EC1(I,II-1))  GO  TO  42 

END  IF 
:      IF  (II.GE.7)  GO  TO  42 
DO  23  1=1, NIT 

TN(I)  =  DM(I)  *  R(I)'lc*Cia/II)  /  (  S1(I. II)  -  C1(I,II)  +  R(I) 

&       *d  -  R(i)**ci(i,iih  /  (i  -R(i)5) 

23  CONTINUE 
SUM=0 . 
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DO  24  1=1, NIT 

24  SUH=SU>!+TN(I) 
DO  25  1=1. NIT 

25  ua)=SUH-TN(I) 
DC  25  1=1, MIT 

26  R(I)=DM(I)/(DM(I)+U(I)) 

11=11+1 
SEC1=0. 
GO  TO  41 
C 
42   WRITE(7,503) 
DO  27  1=1, NIT 
C(I)=C1(I,II' 
S(I)=S1(I,II 


IPC(I)=INT(C(I)+.5) 
IPS(I)=INT(S(I)+.5) 


27  WRITE(7, 504)1, IPC(I),IPS(I) 
WRITE(7/200) 

C 

C     *  COMPUTE  REORDER  POINT: JOINT  ORDER  POLICY 

C 

DO  28  1=1, NIT 

R(I)=RHO(I.II) 
IF(KEY.EQ.1)THEN 

CALL  MUST1(R(I),C(I),DLT(I),SVL,MBUY) 
MR(I)=MBUY 
C(I)=C(I)+MR(I) 
S(I)=S(I)+MR(I) 
IC(I)=INT(C(I)+.5) 
IS(I  =INT(S(I)+.5) 
EC(I)=EC1(I,II)+(MR(I)-DLT(I))*H(T) 
END  IF 
IF(KEY.E0.2)THEN 

CALL  MUST2(R(I),C(I),S(I),DLT(I),SVL,MBUY) 
MR(I)=MBUY 
C(I)=C(I)+MR(I 
S(I)=S(I)+MR(I 
IC(I)=INT(C(I)+.5' 
IS(l)=INT(S(l)+.5' 

EC(I)=EC1(I,II)+(MR(I)-DLT(I))*H(I) 
END  IF 
WRITE(7,505)I,R(I),MR(I)/IC(I),IS(I),EC(I) 

28  SEC=SEC+EC(I) 
SHAVE=(TVC-SEC)/TVC 
WRITE (7 , 506 ) SEC, SHAVE 
STOP 

C 
100   FORMAT ( IX, /, ' **   SELECTED  OPTIONS  : ' , / , IX, ■ (*  KEY  = ' , 12 , 2X, • *  SVL. 
&  =' ,F4.2,2X, ■*  LEAD  TIME  TYPE  =' ,I2,2X, '  *  DATA= ' ,12, 
&2X, '*  S.METHD=' ,12, IX, ')',//, '  *  KEY=  1  :  SVL.  BY  PROB .  OF  NO  SHOR 
&TAGE  PER  REPLENISHMENT  CYCLE1,/,1         2  :  SVL.  BY  FRACT .  OF  DEM 
A  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF  ',/,'  *  LTT=  1  :  CONSTANT 
&LEADTIME       2  :  RANDOM  (UNIFORM  W/   MEAN)  '/,'  *  DATA=1  :  DATA 
&1    2  :  DATA2     3  :  DATA3     4  :  DATA4  ',/,'  *  SEARCH  METHOD  =1  : 
&  MODIFIED  SILVER  METHOD    2  :  HEURISTIC, 
&///,20X,'***   INDEPENDENT  ORDER  POLICY    *** ' 

&,//,lX, 'ITEM' ,3X, 'LT  DEMAND(LT.) ' ,2X, 'EOQ(N-INT.) ' ,2X, ' (INTEGER) ' , 
&2X, 'REORDER1 ,2X, 'EXP.  TVC($) ' , 2X, ' S/O  COST') 

110  FORMAT(I3,2X,F6.2,2X,F4.2) 

111  FORMAT(I3,2X,F7.3.3X,F5.2,2X,F5.2) 

200   FORMAT (//, 20X, '***   JOINT  ORDER  POLICY    ***',/ ,2X, ' ITEM1 . 5X, ' RHO ' 
& , 8X , ' REORDER ' , 5X , ' CANBUY ' , 5X , ■ ORDERUP ■ , 5X , ' EXPECTED  COST ' ) 

501  F0RMAT(1X,I3,4X,F7.3, ' ( ' ,F5 .3 , ' ) ' ,3X, F8.3 , 6X, 14, 7X, 13 ,4X, F9 .3 ,3X, 
83) 

502  FORMAT (3 IX, 'SUM  OF  AVG.  ANNUAL  TVC($)= ' . F9 . 3) 

503  FORMAT( ' ' ,//,15X, ■ **   OPTIMAL  PAIR  (C,S)   *^  ',/, 
&16X, 'ITEM' ,7X, 'OPT.  C',7X.'OPT.  S1) 

504  FORMAT(16X,I3,9X,I4,9X,I4) 

505  FORMAT(2X,I3,4X,F6.4,7X,I4,9X,I4,7X,I4,7X,F10.4) 
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c 
c 


506  FORMAT ( 4 5X, 'TOTAL  COST  =' ,F10 .4 , ' ( ' , F6 .4 , ' )  '  ) 

507  FORMAT ('  '  ,//,3X,  »*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
&,/,15X, ' (WHEN  LEADTIME  =  0) ' ,/ , 15X, '*  AT  ITERATION  :',I2,/, 

&1X , ' ITEM '  .4X  ' RHO ■ , 4X , ' OPT . C ' , 2X , ' OPT . S ' , 3X , ' EXP . COST ' ) 

508  FORMAT ( 2X , I 2 , 3X , F6 . 4 , 4X , 13 , 4X , 13 , 4X , F8 . 4 ) 

509  FORMAT ( 2 3X, 'TOTAL  = ' ,F9 .4, ' ( ' ,F6 .4, ' ) ' ) 
END 

SUBROUTINE  POISON  (KEY,DLT/Q,SVL,MX,M/MAX,CCDF1) 

PSO=l-SVL 
IF(KEY.EQ.l)  THEN 

CALL  SERV1(DLT,PS0,MAX,M) 
RETURN 
END  IF 
IF (KEY. EQ. 2)  THEN 

CALL  SERV2(Q,DLT,PS0,MAX,M,CCDF1) 
RETURN 
END  IF 
END 


SUBROUTINE  MAXD(DLT ,MAX) 

DIMENSION  PMF(750) 

1=1 

PMF(1)=EXP(-DLT) 

CDF=PMF(1) 

GO  TO  2 

1  PMF(I)=(DLT/(I-1))*PMF(I-1) 
CDF=CDF+PMF(I) 

2  IF(CDF.GE. 0.9995)  GO  TO  3 
1=1  +  1 

GO  TO  1 

3  MAX  =  1+1 
RETURN 
END 


C 

c 


SUBROUTINE  SERV1 (DLT ,PSO ,MAX,M) 
DIMENSION  PMF(200),CDF(200),CCDF(200) 


PMF(1)=EXP(-DLT) 
CDF(1)=PMF(1) 
CCDF(1)=1-CDF(1) 
CALL  MAXD(DLT,MAX) 
DO  1  I=2,MAX+1 

PMF(I)=(DLT/(I-1))*PMF(I-1) 

CDF(I)=CDF(I-1)+PMF(I) 

1  CCDF(I)=1-CDF(I) 
DO  2  I=1,MAX+1 

2  IF(PSO.LT.CCDF(I) .AND.PSO.GT. CCDF(I+1) )GO  TO  3 

3  M=I 
RETURN 
END 

SUBROUTINE  SERV2 (Q,DLT,PSO,MAX,M, CCDF1 ) 

DIMENSION  PMF(200),CDF(200),CCDF(200) 

STON=Q^PSO 

PMF(1)=EXP(-DLT) 

CDF(1)=PMF(1) 

CCDF(1)=1.-CDF(1) 

CALL  MAXD(DLT,MAX) 

DO  1  I=2,MAX+1 

PMF(I)=(DLT/(I-1))*PMF(I-1) 

CDF(I)=CDF(I-1)+PMF(I) 

CCDF(I)=1-CDF(I) 

1  CONTINUE 
1=0 

2  IF(I.EQ.O)THEN 
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PST=DLT 
ET  SE 

'^FST=DLT-I-DLT*CDF(I)+I*CDF(It1) 
END  IF 
IF(PST.LT.STON)  THEN 

K=I 

CCDF1=CCDF(M+1) 
ELSE 
1=1+1 
GO  TO  2 
END  IF 
RETURN 
END 

SUBROUTINE  LEAD(TLT,TLY) 


CALL  LRND(719325,U, 1,1,0) 

TLY=2*TLT*U 

RETURN 

END 

SUBROUTINE  MUST1 (R, C,DLT, SVL ,MBUY) 

DIMENSION  PMF(400) ,CDF(400) 

IC=INT(C+.5) 

RC=(1./R)**IC 

RHS=SVL/(R**IC) 

PMF(1)=EXP(-DLT) 

IF(PMF(1) .LT.1.0E-70)  PMF(1)=0. 

CDF(1)=PMF(1) 

DO  1  1=2,400 

PMF(I)=(DLT/(I-1))*PMF(I-1) 
IF(PMF(I) .LT.1.0E-70)  PMF(I)=0. 
CDF(I)=CDF(I-1)+PMF(I) 

1  CONTINUE 
JS=-2 

3     SUM=0 . 
KF=JS+2 
XL=JS+IC 
DO  2  IX=KF,KL 

2  SUM=SUM+PMF(IX+1)*((1/R)**IX) 
TEST=RC*CDF ( KL+1 ) -PMF ( JS+2 ) - (R** ( JS+1 ) ) *SUM 
IF(TEST.GE.RHS)  THEN 

MBUY=JS 
ELSE 

JS=JS+1 

GO  TO  3 
END  IF 
RETURN 
END 

SUBROUTINE  MUST2 (R , C , S , DLT , SVL , MBUY) 

1 

DIMENSION  PMF(400),CDF(400) 

IC=INT(C+.5) 

RHS=(S-C+R*(1-R**C)/(1-R))*(1-SVL) 

PMF(1)=EXP(-DLT) 

IF(PMF(1) .LT.1.0E-70)  PMF(1)=0. 

CDF(1)=PMF(1) 

DO  1  1=2,400 

PMF(I)=(DLT/(I-1))*PMF(I-1) 
IF(PMF(I) .LT.1.0E-70)  PMF(I)=0. 
CDF(I)=CDF(I-1)+PMF(I) 
1     CONTINUE 

JS=-1 

3  SUM=0 . 

DO  2  IW=JS+1,JS+IC 
IF(IW.LT.l)  THEN 
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c 
c 


VCDF=0 . 

IF(  IW.EQ.0)VCDF1=CDF(1) 
ELSE 
VCDF=CDF(IW) 
VCDF1=CDF(IW+1) 
END  IF 
SUM  =  SUM+(1/R**IW)*(DLT-IW-DLT*VCDF+IW*VCDF1) 
IF(JS.LT.l)  THEN 
VCDF=0. 

IF(  JS.EQ.0)VCDF1=CDF(1) 
ELSE 

VCDF=CDF(JS) 
VCDF1=GDF(JS+1) 
END  IF 

TEST=R**C*(DLT-JS-DLT*VCDF+JS*VCDF1)+(1-R)*R**(C+JS)*SUM 
IF(TEST.LE.RHS)  THEN 

MBUY=JS 
ELSE 

JS=JS+1 
GO  TO  3 
END  IF 
RETURN 
END 

SUBROUTINE  OPTION  (ID ,LT ,KEY,LTT ,SVL ,MSRCH) 


WRITE(6,50) 
READ(5,51)  ID 
WRITE (6, 53) 
READ (5, 51)  KEY 
WRITE(6,55) 
READ (5, 57)  SVL 
WRITE(6,52) 
READ(5,51)  LT 
WRITE(6,54) 
READ (5, 51)  LTT 
WRITE (6, 56) 
RE AD ^5, 51)  MSRCH 
RETURN  ' 

'50   FORMATCl',1  *  ENTER  OPTION  FOR  DATA  FILE  (  1-  4)  !  ',/,'     1  = 
&DATA1  2  =  DATA2   3  =DATA3   4  =  DATA4 • ) 

51  FORMAT (II) 

52  FORMAT ('l',1  *  ENTER  THE  LEADTIME  IN  MONTH  !  ',/,5X,  ' 1  =   1  MONTH 
&',/,5X,'3  ■   3  MONTH' ,/,5X, '6  =   6  MONTH1) 

53  FORMAT ('1','  *  ENTER  OPTION  FOR  THE  METHOD  OF  SVL.  MEASURE  !  ',/,' 
&     1  =   PROB.  OF  NO  SHORTAGE  PER  REPLENISHMENT  CYCLE  \/,4X,'  2 
&=   FRACTION  OF  DEMAND  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF  ' ) 

54  FORMAT ( ' 1 ' , '  *  ENTER  THE  LEADTIME  TYPE     !  ' ,/,5X,  ' 1  =  CONSTANT 
&  ' ,/,5X,'2  =  RANDOM  (UNIFORM  W/  MEAN)') 

55  FORMATCl1,'  *  ENTER  THE  DESIRED  SERVICE  LEVEL   !  ',/. 4X, ' (SHOULD 
&BE  IN  THE  FORM  OF  "F4.2"  :  EX.  ;  0.90  OR  0.99  ...ETC.)') 

56  FORMAT ('1','  *  ENTER  OPTION  FOR  SEARCH  METHOD  (1-2 )!',/,'     1  = 
&MODIFIED  SILVER  METHOD    2  =  HEURISTIC) 

57  FORMAT (F4. 2) 
END 


INPUT  FORMAT 


004 

050.00   C 

).20 

1 

290.000 

6, 

,90 

10, 

.00 

2 

41.000 

1, 

.20 

10, 

.00 

3 

77.000 

■3 

.90 

10, 

.00 

4 

122.000 

2. 

.30 

10, 

.00 
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2.  EXEC  PROGRAM 

&TRACE  OFF 

&FN  =  PARA 

&FNO  =  &CONCAT  OF  &FN  OUTPUT 

&TYPE  Do  you  need  to  compile  your  program  ?  (Y) 

&READ  VAR  &R_COMPILE 

&IF  &R_COMPILE  NE  Y  &GOTO  -RUN 

-H  FORTVS  &FN 

&IF  &RC  EQ  0  &GOTO  -RUN 

&TYPE  Your  program  did  not  compile;  check  for  errors. 

&TYPE  Do  you  wish  to  view  the  program  LISTING  file?  (Y) 

&READ  VAR  &RSP1 

&IF  &RSP1  EQ  Y  BROWSE  &FN  LISTING  A 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y) 

&READ  VAR  &RESP1 

&IF  &RESP1  NE  Y  &EXIT  1 

&COMMAND  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP2 

&IF  &RESP2  EQ  Y  &GOTO  -H 

&EXIT  1 

-RUN 

FILEDEF  01  DISK  &FN  DATA1  Al 

FILEDEF  02  DISK  &FN  DATA2  Al 

FILEDEF  03  DISK  &FN  DATA3  Al 

FILEDEF  04  DISK  &FN  DATA4  Al 

FILEDEF  07  DISK  &FN  OUTPUT  Al  (LRECL  133 

LOAD  &FN  (START 

&IF  &RC  EQ  0  &SKIP  9 

&TYPE  Your  program  did  not  run  correctly;  check  for  errors 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y) 

&READ  VAR  &RESP3 

&IF  &RESP3  NE  Y  &EXIT  2 

&C0MMAND  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP4 

&IF  &RESP4  EQ  Y  &GOTO  -H 

&EXIT  2 

&TYPE  YOUR  OUTPUT  IS  IN  THE  FILE  &FN  OUTPUT  A 

&TYPE  Do  you  wish  to  BROWSE  your  output?  (Y) 

&READ  VAR  &RESP 

&IF  &RESP  EQ  Y  &C0MMAND  3R0WSE  &FN  OUTPUT  A 

&TYPE  Print  your  output  file?   (Y) 

&READ  VAR  &RESP7 

&IF  &RESP7  EQ  Y  &COMMAND  PRINT  &FN  OUTPUT  A 

-REDO 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y/N) 

&READ  VAR  &RESP5 

&IF  &RESP5  EQ  Y  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP6 

&RESP56  =  &CONCAT  OF  &RESP5  &RESP6 

&IF  &RESP56  EQ  YY  &GOTO  -H 

&IF  &RESP6  EQ  Y  &GOTO  -RUN 

&EXIT 
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3.  RUN  EXAMPLE 

Sara 
o  you  need  to  compile  your  program  ?  (Y) 

VS  FORTRAN  COMPILER  ENTERED.   11:58:02 

**MAIN**  END  OF  COMPILATION  1  ****** 

**POISON**  END  OF  COMPILATION  2  ****** 

**MAXD**  END  OF  COMPILATION  3  ****** 

**SERV1**  END  OF  COMPILATION  4  ****** 

**SERV2**  END  OF  COMPILATION  5  ****** 

**LEAD**  END  OF  COMPILATION  6  ****** 

**MUST1**  END  OF  COMPILATION  7  ****** 

**MUST2**  END  OF  COMPILATION  8  ****** 

**OPTION**  END  OF  COMPILATION  9  ****** 
VS  FORTRAN  COMPILER  EXITED.    11:58:10 

EXECUTION  BEGINS... 

*  ENTER  OPTION  FOR  DATA  FILE  (1-4)  ! 

1  =   DATA1  2  =  DATA2   3  =DATA3   4  =  DATA4 
1 

*  ENTER  THE  LEADTIME  IN  MONTH  ! 

1=1  MONTH 
3=3  MONTH 
6=6  MONTH 
1 

*  ENTER  OPTION  FOR  THE  METHOD  OF  SVL.  MEASURE  ! 

1  =   PROB.  OF  NO  SHORTAGE  PER  REPLENISHMENT  CYCLE 

2  =   FRACTION  OF  DEMAND  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF 
2 

*  ENTER  THE  LEADTIME  TYPE     ! 

1  =  CONSTANT 

2  =  RANDOM  (UNIFORM  W/  MEAN) 
1 

*  ENTER  THE  DESIRED  SERVICE  LEVEL   ! 

(SHOULD  BE  IN  THE  FORM  OF  "F4.2"  :  EX.  ;  0.90  OR  0.99  ...ETC.) 
0.99 

*  ENTER  OPTION  FOR  SEARCH  METHOD  (1-2)! 

1  =   MODIFIED  SILVER  METHOD    2  =  HEURISTIC 
1 

YOUR  OUTPUT  IS  IN  THE  FILE  PARA  OUTPUT  A 
Do  you  wish  to  BROWSE  your  output?  (Y) 


rint  your  output  file?  (Y) 
n 
Do  you  wish  to  XEDIT  the  program  file?  (Y/N) 

n 

Do  you  wish  to  run  the  program  again?  (Y) 

R;  T=1.47/2.39  11:58:39 
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4.  OUTPUT  1  (  P2  =  0.93  ) 


**   SELECTED  OPTIONS  : 
(*  KEY  =2   *  SVL.  =0.99 


*  LEAD  TIME  TYPE  =  1   *  DATA=  1   *  S.METHD=  1  ) 


*  KEY=  1  :  SVL.  BY  PROB.  OF  NO  SHORTAGE  PER  REPLENISHMENT  CYCLE 

2  :  SVL.  BY  FRACT.  OF  DEMAND  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF 

*  LTT=  1  :  CONSTANT  LEADTIME       2  :  RANDOM  (UNIFORM  W/   MEAN) 

*  DATA=1  :  DATA1    2  :  DATA2     3  :  DAT A3     4  :  DATA4 

*  SEARCH  METHOD  =1  :  MODIFIED  SILVER  METHOD    2  :  HEURISTIC 

***   INDEPENDENT  ORDER  POLICY    *** 


ITEM    LT  DEMAND (LT.] 

1  24.167(0.083 

2  3.417(0.083 

3  6.417(0.083 

4  10.167(0.083) 

i   EOQ(N-INT.)   (INTEGER)   REORDER 
i    158.800       159        25 
I     143.178       143         3 
i     108.340       109         7 
i     178.399       178        10 

SUM  OF  AVG.  ANNUAL  TVC($)  = 

EXP.  TVC($) 

220.984 

34.383 

35.740 

32.217 

=   423.323 

S/O 

1. 
1, 

3 

1 , 

COST 
.983 
.382 
.499 
.536 

*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  =  0) 
*  AT  ITERATION  :  1 
ITEM    RHO     OPT.C   OPT . 5    EXP . COST 

1  0.9942      79     146     201.3498 

2  0.9272     31     90     21.7035 

3  0.9650     27     85     66.6963 

4  0.9774     45     138     63.6090 

• 

TOTAL  =  353.3584(0.1641) 


*  PARAMETERS  &   EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  =  0) 


*  AT  ITERATION 

:  2 

ITEM 

RHO 

OPT.C   OPT.S 

EXP. COST 

1 

0.9972 

84    152 

210.4767 

2 

0.9499 

36     97 

23.4369 

3 

0.9768 

48     89 

70.0377 

4 

0.9850 

71    145 

66.8895 

TOTAL  = 

370.3406(0 

1240) 


*  PARAMETERS  5c  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  =  0) 


*  AT  ITERATION 

:  3 

ITEM 

RHO 

OPT.C   OPT.S 

EXP. COST 

1 

0.9973 

34     152 

210.5555 

2 

0.9464 

36     96 

23.1011 

3 

0.9741 

43     88 

68.8164 

4 

0.9835 

70    143 

65.8673 

TOTAL  =  368.3406(0.1299) 


*  PARAMETERS  &   EX. COSTS  FOR  JOINT  REPLENISHMENT 


(WHEN  LEADTIME 
*  AT  ITERATION 

=  0) 

:  4 

ITEM 

RHO 

OPT.C   OPT.S 

EXP. COST 

1 

0, 

.9974 

87     152 

210.8846 

2 

0 

.9469 

36     96 

23.1429 

3 

0 

.9744 

43     88 

68.9584 

4 

0 

.9835 

70     143 

65.9083 

TOTAL  =  363.3940(0.1286 
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* 

PARAMETERS 

£  EX. COSTS  FOR 

JOINT  REPLENISH^ 

(WHEN  LEADTIME  ' 

=  0) 

*  AT  ITERATION 

:  5 

El 

!    RHO 

CFT.C   OPT.S 

EXP. COST 

1 

J. 

0.9973 

87     152 

210.8239 

2 

0.9467 

36     96 

23.1271 

3 

0.9744 

43     88 

68.9309 

4 

0.9835 

70    143 

65.8689 

TOTAL  = 

368.7505(0.1289) 

**   OPTIMAL  PAIR  (C,S)  ** 

ITEM       OPT.  C  OPT.  S 

1  37  152 

2  36  96 

3  43  88 

4  70  143 


EM 

RHO 

1 

0.9973 

2 

0.9467 

3 

0.9744 

4 

0.9835 

***   JOINT  ORDER  POLICY  *** 

REORDER     CANBUY  ORDERUP 

25          112  177 

-1          35  95 

5          48  93 

7           77  15Q 


EXPECTED  COST 

211.9739 

22.0671 

67.8259 

64.4122 


TOTAL  COST  =   366.2788(0.1348) 
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5.  OUTPUT  2  (PI  =  0.95  ) 


**   SELECTED  OPTIONS  : 
(*  KEY  =1   *  SVL.  =0.95 


*  LEAD  TIME  TYPE  =  1   *  DATA=  1   *  S.METHD=  1  ) 


*  KEY=  1  :  SVL.  BY  PROB .  OF  NO  SHORTAGE  PER  REPLENISHMENT  CYCLE 

2  :  SVL.  BY  FRACT.  OF  DEMAND  TO  BE  SATISFIED  DIRECTLY  FROM  SHELF 

*  LTT=  1  :  CONSTANT  LEADTIME       2  :  RANDOM  (UNIFORM  W/   MEAN) 

*  DATA=1  :  DATA1    2  :  DATA2     3  :  DATA3     4  :  DATA4 


*  SEARCH  METHOD  =1  :  MODIFIED  SILVER  METHOD 


HEURISTIC 


***   INDEPENDENT  ORDER  POLICY   *** 


ITEM 

LT  DEMAND (LT' 

I   EOQ(N-INT.) 

1 

24.167(0.083 

1     158.800 

2 

3.417(0.083 

)     143.178 

3 

6.417(0.083 

I    108.840 

4 

10.167(0.083 

\     178.399 

INTEGER) 

REORDER 

EXP.  TVC($) 

S/O  COST 

159 

33 

232.024 

15.113 

143 

7 

35.343 

16.762 

109 

11 

38.860 

22.051 

178 

16 

84.977 

13.453 

SUM  OF  AVG.  ANNUAL  TVC($)=   441.203 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  =  0) 


*  AT  ITERATION 

:  1 

ITEM 

RHO 

OPT.C   OPT.S 

EXP, 

.COST 

1 

0.9942 

79     146 

201, 

.3498 

2 

0.9272 

31      90 

21 

.7035 

3 

0.9650 

27     35 

66 

.6963 

4 

0.9774 

45     138 

63 

.6090 

TOTAL  =  353.3584(0.1980) 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 


(WHEN  LEADTIME  =  0) 
*  AT  ITERATION 


ITEM 

RHO 

OPT.C 

OPT.S 

EXP. COST 

1 

0.9972 

84 

152 

210.4767 

2 

0.9499 

36 

97 

23.4369 

3 

0.9768 

48 

89 

70.0377 

4 

0.9850 

71 

145 

66.3895 

TOTAL  =  370.3406(0.1595) 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 


(WHEN  LEADTIME  =  0) 


*  AT  ITERATION 

:  3 

ITEM 

RHO 

OPT.C   OPT.S 

EXP. COST 

1 

0.9973 

84    152 

210.5555 

2 

0.9464 

36     96 

23.1011 

3 

0.9741 

43     88 

68.8164 

4 

0.9835 

70    143 

65.8678 

TOTAL  =  368.3406(0.1651) 


*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 
(WHEN  LEADTIME  =  0) 
*  AT  ITERATION  :  4 


ITEM 

RHO 

OPT.C 

OPT.S 

EXP. COST 

1 

0.9974 

87 

152 

210.8846 

2 

0.9469 

36 

96 

23.1429 

3 

0.9744 

43 

88 

68.9584 

4 

0.9835 

70 

143 

65.9083 

TOTAL  =  363.3940(0.1639) 
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*  PARAMETERS  &  EX. COSTS  FOR  JOINT  REPLENISHMENT 


(WHEN  LEADTIME 
*  AT  ITERATION 

=  0) 

:  5 

ITEM 

RHO 

OPT.C  OPT.S 

EXP. COST 

1 

0.9973 

87     152 

210.8239 

2 

0.9467 

36     96 

23.1271 

3 

0.9744 

43     88 

68.9309 

4 

0.9835 

70    143 

65.8689 

TOTAL  =  368.7505(0.1642) 


**   OPTIMAL  PAIR  (C,S)  ** 

ITEM       OPT.  C  OPT.  S 

1  87  152 

2  36  96 

3  43  88 

4  70  143 


***   JOINT 

ORDER  POLICY 

*** 

ITEM 

RHO 

REORDER 

CANBUY 

ORDERUP 

EXPECTED  COST 

1 

0.9973 

32 

119 

184 

221.6339 

2 

0.9467 

4 

40 

100 

23.2671 

3 

0.9744 

9 

52 

97 

70.9459 

4 

0.9835 

13 

83 

156 

67.1722 

TOTAL  COST  =   383.0190(0.1319) 
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APPENDIX  B 
"SUBPROGRAM  TO  TEST  THE  MODEL 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 


VARIABLE  DEFINITIONS 


KEY   =  OPTIONS  FOR  ORDER  POLICY 

ISTA  =  OPTIONS  FOR  START  INVENTORY  LEVEL 

IPRINT=  OPTIONS  FOR  PRINTOUT 

NIT   =  NO. OF  ITEMS  IN  THE  GROUP  TO  BE  EVALUATED 

NSIM  =  SIMULATION  PERIOD 

SIMT  =  SIMULATION  INTERVAL  TIME  (1YEAR) 

DM   =  ANNUAL  DEMAND  RATE 

ME    =  ORDER-UP  POINT 

IEQ   =  ECONOMIC  ORDER  QUANTITY ( INTEGER) 

MC    =  CAN  ORDER  POINT  (CAN  3UY) 

MR   =  REORDER  POINT(Ml=IND. :M2=JOINT) 

TDA   =  TIME  OF  DEMAND  OCCURING  TIME 

TOA   =  ORDER  QTY.  ARRIVAL  TIME 

M    =  INDEX  OF  EVENT  FOR  DEMAND  OCCURED  ITEM 

IA    =  INDEX  OF  EVENT  FOR  OREDR  QTY.  ARRIVED  ITEM 

IET   =  INDEX  OF  EVENT  TYPE 

CT    =  MASTER  CLOCK  TIME  IN  YRS(THE  EARLIEST  EVENT  OCCURING  TIME) 

DA    =  TIME  INTERVAL  BETWEEN  2  SUCCESSIVE  DEMANDS 

ARV   =  UNIFORM  RANDOM  VARIABLE  FOR  LEADTIME  WITH  MEAN   "ULT" 

ULT   =  EXPECTED  LEADTIME 

IQ    =  ORDER  QTY. 

IT    =  INDEX  FOR  REPLENISH  TRIGGER  ITEM 

JT    =  INDEX  FOR  JOINT  REPLENISHED  ITEM 

IOH   =  ONHAND  INVENTORY 

IOHS  =  START  LEVEL 

UOH   =  UNIT  YEARS  ONHAND  QTY 

UBO   =  UNIT  YEARS  BACKORDER  QTY 

NBQ   =  NO.  OF  BACKORDERS  (QTY.) 

NTBQ  =  TOTAL  NUMBER  OF  BACKORDERS 

BCT   =  BACKORDER  TIME 

TBOT  =  TOTAL  BACKORDER  TIME 

NOD   =  TOATL  NUMBER  OF  ORDERS  FOR  EACH  ITEM 

NDM   =  TOTAL  NUMBER  OF  DEMANDS  FOR  EACH  ITEM 

NTDM  =  TOTAL  NUMBER  OF  DEMANDS  DURING  TMAX 

NJTO  =  NUMBER  OF  JOINT  REPLENISHMENTS 

UP  =  UNIT  PRICE 

F  =  HOLDING  COST  RATE 

HC  =  HOLDING  COST 

OC  =  ORDER  &  REVIEW   COST 

OCI  =  INDIVIDUAL  ORDERING  COST 

OCG  =  JOINT  (GROUP)  ORDERING  COST 

SCU  =  BACKORDER  COST  PER  BACKORDER  UNIT 

SCT  =  BACKORDER  COST  (TIME  DEPENDENT  FACTOR) 

TVC  =  AVG.  ANNUAL  TOTAL  VARIABLE  COSTS 


DIMENSION  ISEED(IO) ,DM(10) ,MR( 10 ) ,MC( 10) ,ME ( 10) , IOH( 10) ,ASVL(10) , 
&TDA(100) ,IT(10,100).JT(10,100)/IQ(100),DA(100)/TOA(100),IEQ(10)/ 
&H(10),OCI(10) ,UP(10).SCU(10),SCT(10),F(10) .NJTO(IOO) , OCT (100) , 
&UOH(10,100) ,AUOH(100),UBO(10,100) , AUBO( 100 ) ,NOD( 10 , 100 ) ,KB(10,100) 
&,ANOD(100) ,NDM(10,100) ,ANDM(100) ,NBQ(10,100) ,NTB0( 10 , 100 ) ,KBT(10) , 
&BOT(10,100) , TBOT (10. 100) ,FRC(10,100) , FRO( 10 , 100 ) 7SVL( 10 , 100 ) , 
SAFRQ(IOO) ,OC(10,100) ,AOC(100) , COST( 10 , 100 ) ,TCOST( 100 ) ,NTOD(100) , 
&HCOST(10,100) ,SCOST(10/100),SCTU(10,100),SCTT(10/100) ,NTDM(100), 
&AHCOST(100) ,ASCTU(100) ,ASCTT(100) ,ABOT(10) ,ABOQ(100) ,ACOST(100) , 
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&THCOST(100) ,TSCTU(100) .TSCTT(iOO) ,OCTT(100) ,M1(10) ,M2(1G) , 
&TUCH(100)  ,TUBOUOO)  ,T?RTUO)  ,TFRQ(100) ,NTDMS(100)  ,NT0DS(1C0) , 
SANS (10) 
COMMON  ICASE^SEED^I^DA^RV^LTM^STP^I^ME^EY^EQ 

C 

C   DATA  INITIALIZATION  

C 


c 


SIMT 

= 

1. 

TMAX 

s 

0. 

IC 

= 

0 

ISD 

= 

293715 

CT 

= 

0. 

DT 

= 

0. 

AT 

= 

0. 

K 

= 

1 

N 

= 

0 

TVC 

= 

0. 

NTJTO 

= 

0 

TANDM 

s 

0. 

TANOD 

= 

0. 

TAOC 

= 

0. 

TCST 

= 

0. 

STV 

= 

0. 

ANJ 

= 

0. 

c 

C READ  PROCESS.  OPTIONS 

C 

CALL  OPTION ( ID , RATE1 , RATE2 , KEY , LTOPT , NSIM , ISTA , IP ) 
C 
C READ  INPUT  DATA   

C 

READ ( ID , 101 )NIT , OCG , ISV , MSV , LT 
ULT=LT/12. 
DO  401  1=1, NIT 

READ(ID,100)DM(I),M1(I),M2(I),MC(I),ME(I),IEQ(I),OCI(I),UP(I),F(I) 
IF(KEY.E0.1)MR(I)=Ml(r 
IF(KEY.E0.2)MR(I)=M2(I 
SCU(I)=RATE1*UP(I' 
SCT(I)=RATE2*UP(I 
H(I)=UP(I)*F(I) 
401   CONTINUE 


WRITE (7, 198) ID, KEY, NSIM, ISTA, IP, LTOPT, MSV, ISV, ULT,LT 

DO   1  1=1, NIT 

ISEED(I)=ISD+I**3 

KBT(I)   =  0 

ANB(I)   =0. 

TFRT(I)  =  0. 

ABOT(I   =  0. 

ASVL(I)  =  0. 

TOA(I)   =  10.**10 

TDA(I)   =  0. 
DO   2  1=1, NIT 
DO   2  J=l,100 

DA(J)       =  0. 

IQ(J)        =  0 

SVL(I,J)     =0. 

UOH(I,J)     =  0. 

TUOH(J)      =  0. 

AUOH(J)      =  0. 

COST(I,J)    =  0. 

TCOST(J)     =  0. 

ACOST(J)     =  0. 

UBO(I,J)     =  0. 

TUBO(J      =  0. 

AUBO(J)      =  0. 

NOD(I/J)     =  0 
NTOD(J)      =  0 
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NTODS(J) 

at 

0 

ANOD(J) 

a 

0. 

NDM(I,J) 

s 

0 

NTDM(J) 

= 

0 

NTDMS(J) 

s 

0 

ANDM(J) 

= 

0. 

NBQ(I,J) 

s 

0 

NTSQ(I,J) 
BOTTI.J) 

s 

0 

= 

0. 

TBOT(I.J) 

= 

0. 

ABOQ(J) 

= 

0. 

FRQ(I,J) 
FRC(I,J) 

= 

0. 

= 

0. 

TFRQ(J) 

= 

0. 

AFRO (J) 

= 

0. 

OC(I,J) 

= 

0. 

OCT(J) 

= 

0. 

OCTT(J) 

= 

0. 

AOC(J) 

= 

0. 

NJTO(J) 

s 

0 

HCOST(I,J) 

= 

0. 

THCOST(J) 

= 

0. 

AHCOST(J) 

= 

0. 

SCOST(I,J) 

= 

0. 

SCTU(I,J) 

= 

0. 

TSCTU(J) 

= 

0. 

ASCTU(J) 

= 

0. 

SCTT(I,J) 

= 

0. 

TSCTT(J) 

= 

0. 

ASCTT(J) 

= 

0. 

KB   (I, J) 

= 

0 

IT   (I, J) 

= 

0 

2 

JT   (I, J) 

DO  3   1=1, NIT 

= 

0 

IF(ISTA.EQ, 

.l.i 

^ND. KEY. EQ.l) THEN 

ioh(i)  =  : 

CEQC 

END  IF 

IF(ISTA.EQ 

,l.i 

^ND.KEY.EQ.2)  THEN 

IOH(I)   = 

ME  (I) 

END  IF 

IF(ISTA.EQ, 

•  2) 

THEN 

CALL  START 

IOH(I)=ISTP 

END  IF 

3 
C 
C 

CONTINUE 

PRINTOUT  ] 

[NPUT  DATA   

c 

DO   4  1=1, NIT 

4  WRITE(7,199)  I ,DM(I) , IEQ(I) ,NR(I) ,MC(I) ,ME(I) , IOH(I) ,ULT 
WRITE(7,97) 

DO  5   1=1, NIT 

WRITE(7,98)  I,UP(I),F(I),H(I),OCG,OCI(I)/SCU(I)/SCT(I) 
IC  =  IC  +  1 

IF  (IC.GE.  65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 

WRITE(7,200) 
END  IF 

5  CONTINUE 
C 

C   START  SIMULATION  

C 
C 

IF  (IP.EQ.1.0R.IP.EQ.2)THEN 

WRITE(7,200) 

END  IF 
C 

DO  777  ITER=1,NSIM 

IF  (ITER.GT.1.AND.IP.EQ.1.0R.ITER.GT.1.AND.IP.EQ.2)'  THEN 
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IC  =  0 

WRITE (7 ,200) 
END  IF 
TMAX=TMAX+SIMT 

C 
20    IF(CT.GT.TMAX)  GO  TO  77 
N  =  N+l 
JOINT  =  0 
C 

C  DETERMINE  THE  INITIAL  DEMAND  OCCURANCE  TIMES  FOR  EACH  ITEM  

C 

IF(N.LE.l)  THEN 
DO   6   1=1, NIT 
ICASE=I 
CALL  DEMAND 
TDA(I)=DA(ICASE) 
6  UOH(I,ITER)=UOH(I,ITER)+IOH(I)*DA(ICASE) 

C 

C  DETERMINE  THE  EARLIEST  EVENT   

C     (DEMMAND  OCCURED  /  ORDER  PLACED  /  ORDERED  QTY.  ARRIVED) 


C 


IF(TDA(1).LE.TDA(2))  THEN 

TD=TDA(1) 
ELSE 

TD=TDA(2) 
END  IF 

DO  7  1=1, NIT 
IF(TDA(I).LE.TD)  THEN 

TD=TDA(I) 

M=I 
END  IF 

7  CONTINUE 
C 

IET=0 
CT=TDA(M) 
IOH(M)=IOH(M)-l 

IF(MSV.E0.1.AND.IOH(M) .EQ.-l)  KB(M, ITER)=KB(M, ITER)+1 
NDM(M,ITE"R)=NDM(M,ITER)  +  1 
C 

C   SET  NEXT  DEMAND  OCCURING  TIME  FOR  THE  ITEM  :'M'   

C 

ICASE=M 
CALL  DEMAND 
TDA(M)=CT+DA(ICASE) 
GO  TO  90 
END  IF 
C 

IF(TDA(1).LE.TDA(2))  THEN 

TD=TDA(1) 
ELSE 

TD=TDA(2) 
END  IF 
DO  8  1=1, NIT 

IF(TDA(I).LE.TD)  THEN 
TD=TDA(I) 
M=I 
END  IF 

8  CONTINUE 

IF(TOA(l).LE.TOA(2))  THEN 

TA=TOA(l) 
ELSE 

TA=TOA(2) 
END  IF 
DO  9  1=1, MIT 

IF(TOA(I) .LE.TA)  THEN 
TA=T0A(I) 
IA=I 
END  IF 


C 
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9    CONTINUE 
C 

IF(TA.LE.TD)THEN 
GO  TO  93 

ELSE 
IET=0 
CT=TD 
ICASE=M 
CALL  DEMAND 
TDA(M)=CT+DA(ICASE) 
IOH(M)=IOH(M)-l 

IF(MSV.E0.1.AND.IOH(M) .EQ.-l)  KB(M, ITER)=KB(M, ITER)+1 
NDM(M, ITER )=NDM(M, ITER) +1 
END  IF 
C 

C   DETERMINE  SHOULD  ORDER   

C 
90    IF(IOH(M) .GT.MR(M))THEN 

UOH(M,ITER)=UOH(M,ITER)+IOH(M)*DA(M) 

IF(IP.EQ.1)THEN 

WRITE (7, 501 )CT,IET,M,UOH(M, ITER), (IOH(I), 1=1,4) 

END  IF 

IF(IP.EQ.2)THEN 

WRITE (7, 601 )CT,IET,M,UOH(M, ITER), IOH(M) 

END  IF 

IC  =  IC  +  1 
IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.6-5.AND.IP.EQ.2)  THEN 
IC  =  0 

WRITE(7,200) 
END  IF 
GO  TO  20 
END  IF 
C 

IF(K.EQ.1)THEN 
IET=1 
K=K+1 
IT(M,K)=M 
IF(IOH(M).LT.O)  THEN 
NBQ(M,ITER)=-IOH(M) 

IF(KEY.EQ.l)  IQ(M)  =  IEQ(M)+NBQ(M, ITER) 
IF(KEY.EQ.2)  IQ(M)  =  ME [M)+NBQ(M, ITER) 
IF(LTOPT.EQ.l)TOA(M)  =  CT  +  ULT 
IF(LTOPT.EQ.2)THEN 
CALL  ARRIVE 
TOA(M)=CT+ARV 
END  IF 

BOT(M,ITER)=TOA(M)  -  CT 

IF (NBQ ( M , ITER ) . EO . 1 ) TBOT (M , ITER) =T30T (M , ITER ) +BOT (M , ITER ) 
UBO (M , ITER ) =UBO ( M , ITER ) +BOT (M , ITER ) 
IF(IP.EQ.1)THEN 

WRITE(7, 502)CT,IET,M,NBQ(M, ITER) ,BOT(M, ITER) ,UBO(M, ITER) , 
&  IQ(M),T0A(M),(I0H(I),I=1,4) 

END  IF 
IF(IP.EQ.2)THEN 

WRITE(7,602)CT,IET,M,NBQ(M,ITER) ,BOT(M,ITER) ,UBO(M,ITER) , 
&  IQ(M) ,TOA(M),IOH(M) 

END  IF 

IC  =  IC  +  1 
IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 
WRITE(7,200) 
END  IF 
ELSE 

IF(KEY.EQ.l)  IQ(M)  =  IEQ(M) 
IF(KEY.EQ.2)  IO(M)=ME (M) -IOH(M) 
IF(LTOPT.E0.1)TOA(M)=CT+ULT 
IF(LT0PT.EQ.2)THEN 
CALL  ARRIVE 
TOA(M)=CT+ARV 
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END  IF 
IF(IP.EQ.1)THEN 

WRITE(7,503)CT,IET,M,IQ(M) ,TOA(K) , (IOH(I) ,1=1 ,4) 

END  IF 
IF(IP.EQ.2)THEN 

WRITE(7,603)CT,IET,M,IQ(M) ,TOA(M) , IOH(M) 

END  IF 

IC  =  IC  +  1 

IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 
WRITE (7 ,200) 
END  IF 
END  IF 
NOD(M,ITER)  =  NOD (M, ITER) +1 
OC(M,ITER)  =  CC(M,'lTER)+OCG+OCI(M) 
C 

IF(KEY.EQ.l)GO  TO  20 
IF(KEY.EQ.2)GO  TO  92 
C 

END  IF 
C 

C   DON'T  ORDER  UNTIL  PREVIOUS  ORDER  ARRIVES  --- 

C 

IF(M.EQ.IT(M,K))THEN 
IET=0 

IF(IOH(M).LT.O)  THEN.         .  .  . 

NB0(M,ITER  =-IOH(M) 
BOT (M .ITER ) =TOA (M ) -CT 

IF (NBO (M,  ITER )  .EQ.  l)TBOT(M',  ITER )=TBOT(M,  ITER)  +BOT(M,  ITER) 
UBO  ( M7ITER )  =UBO  ( M ,  ITER )  +BOT  (71 ,  ITER ) 
IF(IP.EQ.1)THEN 

WRITE(7,504)CT,IET,M,NBQ(M,ITER) ,30T(M, ITER) , UBO(M, ITER) , 
&  (I0H(I),I=1,4) 

END  IF 

IF(IP.EQ.2)THEN 

WRITE(7,604)CT,IET,M,NBQ(M,ITER) ,BOT(M,ITER) ,UBO(M, ITER) , IOH(M) 
END  IF 

IC  =  IC  +  1 
IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 
WRITE(7,200) 
END  IF 
ELSE 

UOH(M,ITER)=UOH(M,ITER)+IOH(M)*DA(M) 

IF(IP.EQ.1)THEN 

WRITE(7,501)CT,IET,M,UOH(M,ITER) , (IOH( I ) , 1=1 ,4) 

END  IF 

IF(IP.EQ.2)THEN 

WRITE (7, 601 )CT,IET,H,UOH(M, ITER), IOH(M) 

END  IF 

IC  =  IC  +  1 
IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 
WRITE(7,200) 
END  IF 
END  IF 
GO  TO  20 
ELSE 

DO  10  J=1,K 
10      IF(IQ(IT(M,J)).GT.0.OR.IQ(JT(M,J)).GT.0)GO  TO  92 
IET=1 
K=K+1 
IT(M,K)=M 
IF(IOH(M) .LT.0)THEN 

NBQ(M.ITER)   =  -IOH(M) 

IFtKEY.EO.l)IO(M)  =  IEO^M)+NBQ(M, ITER) 
IF(KEY.EQ.2)l6(M)  =  ME tN)+NBQ(M, ITER) 
IF(LTOPT.EQ.l)TOA(M)  =  CT  +  ULT 
IF(LTOPT.EQ.2)  THEN. 
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CALL  ARRIVE 
TOA(M)=CT+ARV 
END  IF 

BOT(M,ITER)  =  TCA(M)  -  CT 

IF (NBQ(M, ITER). EO . 1 )TBOT(M, ITER)=TEOT(M, ITER)+BOT(M, ITER) 
UBO(M,ITER)=UBO(M\ITER)+BOT(M,ITER) 
IF  (IP.  EQ.l)  THEN- 
WRITE^, 502)CT,IET,M,NBQ(M,  ITER), BOT(M, ITER) ,UBO(M, ITER), 
&  IQ(M),T0A(MJ,(I0H(I),I=1,4) 

END  IF 
IF(IP.EQ.2)THEN 

WRITE (7, 602 )CT, IET,M,NBQ(M, ITER) ,BOT(M, ITER) ,UBO(M,ITER)  , 
&  IQ(M),TOA(M),IOH(M) 

END  IF 

IC  =  IC  +  1 
IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 
WRITE(7,200) 
END  IF 
ELSE 

IF(KEY.EQ.1)IQ(M)=IEQ(M) 
IF(KEY.EQ.2)IQ(M)=ME(M)-IOH(M) 
IF ( LTOPT . EQ . 1 ) TOA(M)=CT+ULT 
IF(LTOPT.EQ.2)THEN 
CALL  ARRIVE 
TOA(M)=CT+ARV 
END  IF 
IF(IP.EQ.1)THEN 

WRITE(7,503)CT,IET,M,IQ(M) ,TOA(M), ( IOH(I ) , 1=1 ,4) 
END  IF 
IF(IP.EQ.2)THEN 

WRITE(7,603)CT,IET,M,IQ(M) ,TOA(M) , IOH(M) 
END  IF 

IC  =  IC  +  1 
IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 
WRITE(7,200) 
END  IF 
END  IF 
NOD(M,ITER)  =  NOD(M,ITER)+l 
OC(M,ITER)  =  OC(M,ITER)+OCG+OCI(M) 
END  IF 

IF(KEY.EQ.l)GO  TO  20 
C 

C   DETERMINE  THE  JOINTLY  REPLENISHABLE  ITEM  

C 
91   DO  11  1=1, NIT 

IF  (IOH(I).GT.MR(I).AND.IOH(I).LE.MC(I).AND.IQ(I).EQ.0)THEN 
JT(I.K)  =  I 
IQ(I)  =  ME(I)-IOH(I) 
IF(LT0PT.EQ.1)T0A(I)=CT+ULT 
IF(LTOPT.EQ.2)THEN 
CALL  ARRIVE 
TOA(I)=CT+ARV 
END  IF 

NODa,ITER)=  NOD(I,ITER)  +  l 
OC(I,ITER)  =  OC(I,ITER)+OCI(I) 
JOINT  =  1 

IF ( IP . EQ . 1 . OR . IP . EQ . 2 ) THEN 
WRITE(7,506)I,IQ(I),TOA(I) 
END  IF 

IC  =  IC  +  1 
IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 
WRITE(7,200) 
END  IF 
END  IF 
11    CONTINUE 

NJTO(ITER)  =  NJTO(ITER)  +JOINT 
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GO  TO  20 

92  IF(IOH(M).LT.O)  THEN 

NBQ(M,ITER)=-IOH(M) 

UBO(M,ITER)=UBO(M,ITER)+BOT(M,ITER) 
BOT(M,ITER)=TOA(M)-CT 

IF (NBQ (M , ITER ) . EQ . 1 ) TBOT (M , ITER ) =TBOT (M , ITER) +BOT (M , ITER) 
IF(IP.EQ.1)THEN 

WRITE(7,504)CT,IET.M,NBQ(M/ITER),BOT(M/ITER)/UBO(M,ITER), 
&  (IOH(I),I=l,4) 

END  IF 

IF(IP.EQ.2)THEN 

WRITE(7,604)CT,IET,M,NBQ(M,ITER) ,BOT(M,ITER) ,UB0(M, ITER) , IOH(M) 
END  IF 

IC  =  IC  +  1 
IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 
WRITE(7,200) 
END  IF 
ELSE 

UOH(M,ITER)=UOH(M,ITER)+IOH(M)*DA(M) 

IF(IP.EQ.1)THEN 

WRITE(7,501)CT,IET,M,UOH(M,ITER)/(IOH(I),I=1,4) 

END  IF 

IF(IP.EQ.2)THEN 

WRITE (7, 601 )CT,IET,H,UOH(M, ITER), IOH(M) 

END  IF 

IC  =  IC  +  1 
IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 
WRITE (7, 200) 
END  IF 
END  IF 
GO  TO  20 
C 

C  SET  NEW  ON  HAND  LEVEL  FOR  THE  ITEM  : ' IA '  

C 

93  CT=TA 
IET=2 

IOH(IA)=IOH(IA)+IQ(IA) 
DO  12  1=1, K 

IF(IA.EQ.IT(IA,I))THEN 
IT(IA.I)=0 
IQ(IA)=0 

NTBQ(IA,ITER)=MTBQ(IA,ITER)+NBQ(IA,ITER) 
BOT(IA,ITER)=0. 
NBQ(IA.ITER)=0 
TOA(IA)=10.**10 
END  IF 
IF(IA.EQ.JT(IA,I))THEN 
JT(IA.K)=0 
IQ(IA)=0 

NTBQ ( I A , ITER) =NTBQ ( I A , ITER ) +NBQ ( IA , ITER) 
BOT(IA,ITER)=0. 
NBQ(IA.ITER)=0 
TOA(IA)  =  10>*10 
END  IF 
12    CONTINUE 

IF(IP.EQ.1)THEN 

WRITE(7,500)CT,IET,IA/(IOH(I)/I=1,4) 
END  IF 
IF(IP.E0.2)THEN 

WRITE(7",600)CT,IET,IA,IOH(IA) 
END  IF 

IC  =  IC  +  1 
IF  (IC.GE.65.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2)  THEN 
IC  =  0 
WRITE(7,200) 
END  IF 
"  GO  TO  20 
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c 

C   COMPUTE  INTER- SUMMARY  (ANNUAL)  

C 

77   :f(ip.ne.4>then 
write (7, 300) iter 

END  IF 

DO  15  1=1. NIT 

NTDM ( ITER ) =NTDM( ITER )+MDM( I, ITER) 

NTOD ( ITER) =NTCD ( ITER) +NOD ( I , ITER) 

OCT ( ITER)=CCT ( ITER)+OC ( I , ITER) 

IF(MSV.EQ.1)THEN 

IF(KB(I,ITER) .GT.0)THEN 

FRC(I,ITER)=KB(I,ITER)*1.0/NOD(I,ITER) 
ELSE 

FRC(I,ITER)=0. 
END  IF 
SVL(I,TTER)=(1.-FRC(I,ITER))*100. 
ELSE 

IF(NTBQ(I,ITER).GT.O)THEN 

FRQ(1,ITER)=NTBQ(I,ITER)*1.0/NDM(I,ITER) 
ELSE 

FRO (I, ITER) =0. 
END  IF 

SVL(I,ITER)=(1.-FRQ(I,ITER))*100. 
END  IF 

15  CONTINUE 

.  IF  (.IP. NE.  4)  THEN  ■  . 

DO  16  1=1, NIT 

WRITE(7, 507)1, NDM(I, ITER) ,NOD(I,ITER) ,UOH(I,ITER) ,UBO(I,ITER) 
&NTBQ(I, ITER) ,TBOT(I, ITER) ,SVL(I, ITER) ,IOH(I) 

16  CONTINUE 
WRITE(7,508)NTDM(ITER)  ,NTOD(ITER) 
WRITE (7, 509 )NJTO( ITER) 

END  IF 
C 

C  COMPUTE  COSTS   

C 

DO  17  1=1, NIT 

HCOST ( I , ITER ) =UOH ( I , ITER ) *H ( I ) 

SCTT(I,ITER)=UBO(I,ITER)*(SCT(I)+H(I)) 

SCTU(I,ITER)=NTBQ(I,ITER)*SCU(l) 

SCOST(I , ITER)=SCTU(I , ITER)+SCTT(I , ITER) 

COST (I , ITER)=HCOST(I , ITER)+SCOST(I , ITER)+OC(I , ITER) 

17  TCOST ( ITER )=TCOST( ITER) +COST( I, ITER) 
IF(IP.NE.4)THEN 

WRITE(7,514) 
DO  18  1=1, NIT 

18  WRITE(7,515)  I ,HCOST(I , ITER) , SCTU(I . ITER) , SCTT(I , ITER) , 
&  OC(I, ITER) ,COST(I, ITER) 

WRITE (7, 516) OCT (ITER) , TCOST (ITER) 

END  IF 

DO  19  1=1, NIT 

19  KBT(I)=KBT(I)+KB(I,ITER) 
777   CONTINUE 

C 

C   COMPUTE  FINAL-RESULTS  (MEAN  AND  STD .  DEVIATIONS)  


C 


WRITE(7,510)NSIM 
DO  21  1=1, NIT 
DO  21  J=1,ITER-1 
NTDMS(I)=NTDMS(I)+NDM(I,J' 
NTODS ( I ) =MT0DS ( I ) +NOD ( I . J ' 
TUOH ( I ) =TUOH ( I ) +UOH ( I , J 


TUBO(l)=TUBO(l)+UBO(l,J 
TFRT ( I ) =TFRT ( I  +TBOT ( I , J ) 


TFRO ( I ) =TFRO ( I ) +NTBO ( I , J ) *1 . 0 
OCTT"  ( I )  =OCTT  ( I )  +OC  ( I ,  J  ) 
THCOST ( I ) =THCOST ( I ) +HCOST ( I , J ) 
TSCTT(I)=TSCTT(I)+SCTT(I,J) 
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21  TSCTU(I)=TSCTU(I)+SCTU(I,J) 

DC  22  ~  =  "  MI^ 

ANDK  ( I )  =NTDKS  ( I )  /TMAX 

ANOD  ( I S  =NTODS  ( I )  /TMAX 

AUOH ( I ) =TUCH ( I ) /TMAX 

AUEO ( I } =TUBO ( I ) /TMAX 

ABOT(I)=TFRT(l)/TMAX  ■ 

ABOQ ( I ) =TFRQ ( I ) /TMAX 

ANB(I)  =KBT(I)*1.0/NTODS(I) 

AFRQ ( I ) =TFRQ ( I ) /NTDMS ( I ) 

IF(MSV.EQ.l)  ASVL(I)=(1.-ANB(I))*100. 

IF(MSV.E0.2)  ASVL(I)=(1--AFRQ(I))*100. 

AOC(I)=OCTT(I)/TMAX 

AHCOST ( I ) =THCOST ( I ) /TMAX 

ASCTU(I)=TSCTU(I)/TMAX 

AS  CTT ( I  =TS  CTT ( I ) / TMAX 

ACOST ( I ) =AHCOST ( I ) +AS  CTU ( I ) +ASCTT ( I ) +AOC ( I ) 

22  TVC=TVOACOST(I) 
DO  23  1=1, ITER- 1 

23  NTJTO=NTJTO+NJTO(I) 
IF(NTJTO.GT.0)ANJ  =  NTJTO/TMAX 
DO  24  I=1,ITER-1 

24  STV=STV+(TCOST(I)-TVC)**2 
STDVC=SORT ( STV/ ( ITER- 2 ) ) 

DO  25  1=1, NIT 

25  WRITE(7,511)I,ANDM(I),ANOD(I),AUOH(I),AUBO(I),ABOQ(I),ABOT(I), 
&  ASVL ( I ) 

DO  26  1=1, NIT 
TANDM=TANDM+ANDM ( I ' 
TANOD=TANOD+ANOD ( I 
TAOC=TAOC+AOC(I) 

26  CONTINUE 

WRI TE ( 7 , 5 1 2 ) TANDM , TANOD 
WRITE(7,513)ANJ 
WRITE(7,517) 
DO  27  1=1, NIT 

27  WRITE(7,515)  I ,AHCOST(I) ,ASCTU(I ) ,ASCTT(I) ,AOC(I) ,ACOST(I) 
WRITE(7,516)TAOC,TVC 

WRITE(7,518)STDVC 
STOP 

97  FORMAT (2X, '  ' ,      ///,'*   INPUT  DATA  FOR  COSTS:1,//,'  ITEM' , 2X, ■ UNIT 
&  PRICE' ,2X, 'H. RATE' ,2X, 'H. COST1 ,3X, '   ORDER  COST ' , 9X, ' B/O  COST', 
&/,36X, 'GROUP' ,6X, ' IND . ' ,4X. 'QTY. ' ,5X, 'TIME' ) 

98  FORMAT(2X,I2,4X,7(F6.2,3X)) 

100  FORMAT (F7. 3, 2X, 5(13, 2X) ,3(F6.2,2X)) 

101  FORMAT(I3,3X,F6.2,3X,I2,3X,2(Il,3X)) 

198  FORMAT (//, 5X, ' **  INPUT  DATA  : ' ,/ ,3X,/ .12X, •*  USED  DATA  FILE  NO.  =' 
&,I3,/,12X, '*  ORDER  POLICY  OPTION  = ' , 13 , / , 12X, ' *  STIMULATION  LIMIT 
&(YRS)  =',13  ,  /,12Xf'*  START  LEVEL  OPTION  = ' , 13 , / ,12X, ' *  PRINTOUT 
&OPTION  =' ,I3,/,12X, '*  LEAD  TIME  TYPE   = ' , 13 , / , 12X, /// , 7X, ' *  ORDER 
&POLICY  OPTIONS' ,//,10X, '1  :  INDIVIDUAL  ORDER  POLICY ',/, 10X, ' 2  :  JO 
&INT  ORDER  POLICY' ,//,7X, '*  START  LEVEL  OPTIONS ',//, 10X, ' 1  :  REGULA 
&R  (IND.  POLICY  W/  E.O.Q  ;  JOINT  POLICY  W/   ORDER-UP  QTY. ) ' , / , 10X, ' 
&2  :  RANDOM  (REORDER  POINT  +  1;  ORDER-UP  POINT) ',//, 7X, ' *   PRINTOUT 
&OPTIONS' ,//,10X, '1  :  PRINT  ALL  ITEM ' , / , 10X, ' 2  :  PRINT  ONLY  EVENT 
&ITEM' ,/,10X, '3  :  PRINT  ONLY  SUMMARY ',/, 10X, ' 4  :  PRINT  ONLY  FINAL  R 
&ESULT' ,//,7X, '*  LEAD  TIME  TYPES  \//,10X,'l  :  COSTANT  LEADTIME ' , / , 
&10X,'2  :  RANDOM  LEADTIME  (UNIFORM  W/  MEAN) ' , // , 6X, ' *  SERVICE  LEVEL 
&  :  P',I1,  '=' ,13, IX, '%' ,5X, '*  AVG.  LEAD  TIME  :  ' , F6 . 3 , '  YRS . ' , ' ( ' , 
&I2,'  MON.) ' ,///,3X, 'ITEM' ,3X, 'YR  DEMAND ' ,4X, ' EOQ ', 6X, ' MUST  BUY', 
&3X,'CAN  BUY' ,3X, 'ORD.  UP  QTY' , 3X ,' START  LEVEL'  3X, 'EXP.  LT . ' ) 

199  FORMAT(3X,I2,6X,F7.3.5X,4(I3,3X) , 3X, 13 , 7X, F7 .3) 

200  FORMAT( '1' ,//,35X, • **  SIMULATION  RESULTS  **',//, 20X, ' (EVENT  TYPE 
&  0  :  DEMAND  OCCURED  1  :  ORDER  PLACED  2  :  ORDER  OTY  ARRIVED)' 
&,//,5X, 'EVENT' ,4X, 'ITEM  #' ,2X, 'BACK  ORDER1 ,3X, 'UNIT  YRS (QTY) ' ,4X, ' 
&ON  ORDER ' , 10X , ' ON  HAND ' , 9X , ' JOINT  REPLENISHMENT ' , / , IX, ' TIME ' , 2X , 
&'*TYPE' ,10X, 'QTY. ' ,2X, 'TIME' ,3X, 'ON  HAND' ,2X, 'B/O.S' ,3X, 'QTY. ' ,2X, 
&'AR/T. ' ,4X, '1' ,4X, '2' ,4X, '3' ,4X, '4' ,2X, 'ITEMt' ,2X, 'OD.  QTY. ' ,2X, '0 
&D.AR/T. ' ) 
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c 


c 

c 


300   FORMAT (' 1 ' ,///,20X, '***   INTER  -  SUMMARY   (ANNUAL)   ***',//,  20X, 
&'AT  THE  YEAR  s ' ,2X,I3.// ,3X, • ITEM  #',3X,'#  DEMANDS ' ,3X, ' #  ORDERS' 
&.10X,1  UNIT  YEAR  ' ,7X, ' BACKORDER' ,8X, ' SERVICE ', 6X, ' F/  OH',/,40X, 
&'OM  HAND1 ,3X, '  3/0   ' , 2X; ' QTY. ' ,3X, ' TIME ' , 5X, ' LEVEL(%) ' ) 

500  FORMATS  ■ ,F6.4,2X,X2,3X,X2,48X,4(X4,1X)) 

500  FORMAT ('  ' ,F6 .4,2X, 12 ,3X, 12 , 54X, 14) 

501  FORMATS'  ■ ,FS .4, 2X, 12 ,3X, 12 , 18X, F8 . 2/22X,4(I4, IX) ) 

601  FORMAT{'  ' , F6 .4, 2X, 12, 3X, 12 , 18X,F8 .2 ,28X, 14) 

502  FORMAT('  ' , F6 .4,2X, 12 ,3X, 12 , 5X, 13 ,3X,F6 .3 , 11X,F6 .2 ,3X, 13 , IX, 
&  F6. 3, IX, 4(14, IX)) 

602  FORMAT ('  * , F6 .4, 2X, 12 ,3X, 12 , 5X, 13 ,3X,F6 .3 , 11X,F6 .2 ,3X, 13 , IX, 
&  F6.3,7X,I4) 

503  FORMAT ('  ' , F6 .4, 2X, 12 ,3X, 12 , 37X, 13 , IX, F6 .3 , IX, 4(14, IX) ) 

603  FORMATS  ' , F6 .4 , 2X, 12 , 3X, 12 , 37X , 13 , IX, F6 . 3 , 7X, 14) 

504  FORMAT ('  ' , F6 .4 , 2X, 12 , 3X, 12 , 5X, 13 ,3X, F6 .3 , 11X,F6 .2 , 14X,4(I4 , IX) ) 

604  FORMATS  ' , F6 . 3 , 2X, 12 , 3X, 12 . 5X , 13 , 3X, F6 . 3 , 11X, F6 . 2 , 20X, 14) 

506  FORMAT(86X,I2,6X,I3,5X,F6.3) 

507  FORMAT(5X,I2.9X,I4,7X,I3,9X,F9.4,2X,F6.3,2X,I3,3X,F6.3,5X,F6.2, 
&       9X,I4) 

508  FORMAT (IX, 'TOTAL  = ' , 7X, 15 , 5X, 15 ) 

509  FORMAT (IX, 'TOTAL  NUMBER  OF  JOINT  REPLENISHMENT  =',I5) 

510  FORMAT('l' .///,20X, ' **  FINAL  RESULTS  (AVG.)  :', 10X ,' DURING' , 13 , 2X, 
S'YEARS  ***  '  ,//,3X,  'ITEM  #'  ,3X,  '#  DEMANDS'  ,3X,  '#  ORDERS'  ,9X,  'AVG 
&.  UNIT  YEAR' ,4X, 'AVG.  B/O ', 5X, ' SERVICE '  ,/,40X, 

&'ON  HAND' ,3X, 'B/O.S' ,4X, 'QTY. ' ,2X, 'TIME' ,3X, 'LEVEL' ,'(' ,'%',')' ) 

511  FORMAT(5X,I2,6X,F6.2,5X,F6.2,10X,F7.3,1X,F7.3,1X,F6.2,2X,F7.3, 
-   &2X,F6.2) 

512  FORMAT (IX, 'TOTAL  AVG.= ' , F7 .2 , 5X, F6 .2) 

513  FORMAT (IX, 'AVG.  NUMBER  OF  JOINT  REPLENISHMENT  =',F6.3) 

514  FORMAT(1X,///,10X, '  *  COMPUTED  ANNUAL  COSTS  ($)',//, 2X, 

Sc'ITEM1  ,3X,  'HOLDING'  ,9X,  'BACKORDER'  ,9X,  'ORDER'  ,6X,  'TOTAL'  ,/,22X,  '  QT 
&Y. ' ,7X, 'TIME' ) 

515  FORMAT  2X, 12 , 5X, 5(F7 .3 ,4X) ) 

516  FORMAT ( IX, /,33X, 'TOTAL  =  ' ,2 (F8 . 3 , 3X) ) 

517  FORMAT(1X,///,10X, '  *  COMPUTED  AVG.  ANNUAL  COSTS  ($)',//,2X, 

& ' ITEM ' , 3X , ' HOLDING ' , 9X , ' BACKORDER ' , 9X , ' ORDER ' , 6X , ' TOTAL ' , / , 22X , ' QT 
&Y. ' ,7X, 'TIME' ) 

518  FORMAT ( IX, /,33X, 'STAND.  DEVIATION  =  \F8.3) 
END 

SUBROUTINE  DEMAND 

DIMENSION  ISEED(IO) ,DM(10) ,DA(100) ,MR(10) ,ME(10) , IEQ(IO) 
COMMON  ICASE , ISEED , DM , DA , ARV , ULT , I , ISTP , MR , ME , KEY , IEQ 
CALL  LRND  (ISEED(ICASE) ,U . 1 . 1 , 0 ) 

DA( ICASE )=-(l./DM( ICASE ))*ALOG(U) 
RETURN 
END 

SUBROUTINE  ARRIVE 

DIMENSION  ISEED(IO) ,DM(10) ,DA(100) ,MR(10) ,ME(10) , IEQ(IO) 
COMMON  ICASE , ISEED , DM , DA , ARV , ULT , I , ISTP , MR , ME , KEY , IEQ 
CALL  LRND(ISEED(I),U, 1,1,0) 
ARV=  2.*ULT*U 

RETURN 
END 

SUBROUTINE  START 

DIMENSION  ISEED(IO) ,DM(10) ,DA(100) ,MR(10) ,ME(10) , IEQ(IO) 
COMMON  ICASE , ISEED , DM , DA , ARV , ULT , I , ISTP , MR , ME , KEY , IEQ 
IF(KEY.EQ.l)  THEN 

CALL  LRND (ISEED ( I ),U. 1.1,0) 

IS=  (IEQ(±)-(MR(I)+1))*U 

ISTP=  IS+MR(I)+1 

RETURN 
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END  IF 

IF (KEY. EQ. 2}  THEN 

CALL  LRND(ISEED(I) ,U, 1.1,0) 
15=  (ME(I)-(MR(I)+i))*U 
ISTP=  IS+MR(I)+1 
RETURN 
END  IF 
END 

SUBROUTINE  OPTION (ID, RATE 1 ,RATE2 ,KEY,LTOPT,NSIM, ISTA, IP) 

WRITE(6,35) 
READ(S,34)  ID 
WRITE (6, 39) 
READ ( 5 , 40 ) RATE 1 , RATE2 
WRITE(6,31) 
READ (5, 34)  KEY 
WRITE (6, 37) 
READ (5, 34)  LTOPT 
WRITE(6,36) 
READ (5, 30)  NSIM 
WRITE(6,32) 
READ (5, 34)  ISTA 
WRITE(6,33) 
READ (5, 34)  IP 
RETURN 

30  FORMAT (12) 

31  FORMAT ('l1,1  *  ENTER  OPTION  FOR  ORDER  POLICY  (1  -  2)  !  ' ,/,5X,  ' 
&1  =  INDIVIDUAL  ORDER  POLICY1,/,  5X,'2  =  JOINT  ORDER  POLICY  '  ) 

32  FORMAT ('l1,'   *  ENTER  OPTION   FOR  START  LEVEL  (1  -  2)  !  \/,5X,'l 
&=  REGULAR  (IND. ORDER  POLICY  :  E.O.Q  /  JOINT  ORDER  POLICY  :  ORDER  U 
&P  QTY.) ' ,/.5X, '2  =  RANDOM  BETWEEN  (MUST  ORDER  LEVEL+1 )  AND  (ORDER 
&UP  LEVEL)' ) 

33  FORMAT ('l',1  *  ENTER  OPTION  FOR  PRINTOUT (1  -  4)  !  ■ , / , 5X, ■ 1  =  PRI 
&NT  ALL  ITEM' ,/,5X, '2  =  PRINT  ONLY  EVENT  ITEM'  /,5X,'3  =  PRINT  ONL 
&Y  SUMMARY  \/,5X,'4  =  PRINT  ONLY  FINAL  RESULT') 

34  FORMAT (II) 

35  FORMAT(ll1,1   *  ENTER  OPTION  FOR  DATA  FILE  (1-4)    !  ',/, 
&5X,'l  =  DATA1'  3X,'2  =  DATA2',3X,'3  =  DATA3',3X,,4  =  DATA4 ' ) 

36  FORMAT ('1','  *  ENTER  SIMULATION  LIMIT  YEARS  !',/,5X,'01  =  SIMU 
&LATION  LIMIT  IS  1  YEAR ' , / , 5X, ' 02  =  SIMULATION  LIMIT  IS  2  YEARS' 

&  ,/,7X,'  ',/,5X,'99  =  SIMULATION  LIMIT  IS  99  YEARS  ') 

37  FORMAT('l','  *  ENTER  THE  TYPE  OF  LEADTIME(1  -  2)  !'  /.5X,'l  =  CON 
&STANT  LEADTIME' ,/,5X, '2  =  RANDOM  (UNIFORM  WITH  MEAN)') 

39  FORMATCl1,1   *  ENTER  RATE  OF  STOCKOUT  COST  TO  UNIT  PRICE  !  '  ,  /  ,  5X 
&, '1.25  0.91  =  STOCKOUT  COST  PER  UNIT  IS  125%  OF  UNIT  PRICE ',/, 17X, 
&' STOCKOUT  COST  PER.  TIME  IS  91%  OF  UNIT  PRICE') 

40  F0RMAT(2(F4.2,1X)) 
END 


1.  DATA  FORMAT 
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2.  EXEC  PROGRAM 

&TRACE  OFF 

&FN  =  SIM 

&FNO  =  &CONCAT  OF  ScFN  OUTPUT 

&TYPE  Do  you  need  to  compile  your  program  ?  (Y) 

&READ  VAR  &R_COMPILE 

&IF  &R_COMPILE  NE  Y  &GOTO  -RUN 

-H  FORTVS  &FN 

&IF  &RC  EQ  0  &GOTO  -RUN 

&TYPE  Your  program  did  not  compile;  check  for  errors. 

ScTYPE  Do  you  wish  to  view  the  program  LISTING  file?  (Y) 

ScREAD  VAR  &RSP1 

&IF  &RSP1  EQ  Y  BROWSE  &FN  LISTING  A 

&TYPE  Do  you  wish  to  XEDIT  the  program  file?  (Y) 

ScREAD  VAR  &RESP1 

&IF  &RESP1  NE  Y  &EXIT  1 

ScCOMMAND  XEDIT  &FN  FORTRAN  A 

ScTYPE  Do  you  wish  to  run  the  program  again?  (Y) 

ScREAD  VAR  &RESP2 

ScIF  &RESP2  EQ  Y  &GOTO  -H 

&EXIT  1 

-RUN 

FILEDEF  01  DISK  ScFN  DATA1  Al 

FILEDEF  02  DISK  ScFN  DATA2  Al 

FILEDEF  03  DISK  ScFN  DATA3  Al 

FILEDEF  04  DISK  ScFN  DATA4  Al 

FILEDEF  07  DISK  ScFN  OUTPUT  Al  (LRECL  133 

LOAD  ScFN  (START 

&IF  ScRC  EQ  0  ScSKIP  9 

&TYPE  Your  program  did  not  run  correctly;  check  for  errors 

&TYPE  Do  you  wish  to  XEDIT  the  program  tile?  (Y) 

&READ  VAR  &RESP3 

&IF  &RESP3  NE  Y  &EXIT  2 

&C0MMAND  XEDIT  &FN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP4 

&IF  &RESP4  EQ  Y  &GOTO  -H 

&EXIT  2 

&TYPE  YOUR  OUTPUT  IS  IN  THE  FILE  &FN  OUTPUT  A 

&TYPE  Do  you  wish  to  BROWSE  your  output?  (Y) 

&READ  VAR  &RESP 

&IF  &RESP  EQ  Y  &COMMAND  BROWSE  &FN  OUTPUT  A 

&TYPE  Print  your  output  file?   (Y) 

&READ  VAR  &RESP7 

&IF  &RESP7  EQ  Y  &C0MMAND  PRINT  &FN  OUTPUT  A 

-REDO 

&TYPE  Do  vou  wish  to  XEDIT  the  orogram  file?  (Y/N) 

&READ  VAR"&RESP5 

&IF  &RESP5  EQ  Y  XEDIT  ScFN  FORTRAN  A 

&TYPE  Do  you  wish  to  run  the  program  again?  (Y) 

&READ  VAR  &RESP6 

&RESP56  =  &C0NCAT  OF  &RESP5  &RESP6 

&IF  &RESP56  EQ  YY  &GOTO  -H 

&IF  &RESP6  EQ  Y  &GOTO  -RUN 

&EXIT 
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3.  RUN  EXAMPLE 

sim 

Do  you  need  to  compile  your  program  ?  (Y) 


S  FORTRAN  COMPILER  ENTERED.   12:11:11 

**MAIN**  END  OF  COMPILATION  1  ****** 

**DEMAND**  END  OF  COMPILATION  2  ****** 

**ARRIVE**  END  OF  COMPILATION  3  ****** 

**START**  END  OF  COMPILATION  4  ****** 

**OPTION**  END  OF  COMPILATION  5  ****** 
VS  FORTRAN  COMPILER  EXITED.    12:11:23 

EXECUTION  BEGINS . . . 


f 


*  ENTER  OPTION  FOR  DATA  FILE  (1-4) 

1  =  DATA1    2  =  DATA2    3  =  DATA3   4  =  DATA4 
1 


*  ENTER  RATE  OF  STOCKOUT  COST  TO  UNIT  PRICE  ! 

STOCKOUT  COST  PER  TIME  IS  91%  OF  UNIT  PRICE 


1.25  0.91  =  STOCKOUT  COST  PER  UNIT  IS  125%  OF  UNIT  PRICE 


0.40  0.01 

*  ENTER  OPTION  FOR  ORDER  POLICY  (1  -  2)  ! 

1  =  INDIVIDUAL  ORDER  POLICY 

2  =  JOINT  ORDER  POLICY 


*  ENTER  THE  TYPE  OF  LEADTIME(1  -  2) 

1  =  CONSTANT  LEADTIME 

2  =  RANDOM  (UNIFORM  WITH  MEAN) 


*  ENTER  SIMULATION  LIMIT  YEARS 

01  =  SIMULATION  LIMIT  IS  1  YEAR 

02  =  SIMULATION  LIMIT  IS  2  YEARS 
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99  =  SIMULATION  LIMIT  IS  99  YEARS 

*  ENTER  OPTION   FOR  START  LEVEL  (1  -  2)  ! 

1  =  REGULAR (IND. ORDER  POLICY:  E.0.0  /  JOINT  ORDER  POLICY:  OD.UP  QTY.) 

2  =  RANDOM  BETWEEN  (MUST  ORDER  LEVEL+1)  AND  (ORDER  UP  LEVEL) 
2 

*  ENTER  OPTION  FOR  PRINTOUT (1  -  4)  ! 

1  =  PRINT  ALL  ITEM 

2  =  PRINT  ONLY  EVENT   ITEM 

3  =  PRINT  ONLY  SUMMARY 

4  =  PRINT  ONLY  FINAL  RESULT 
4 

YOUR  OUTPUT  IS  IN  THE  FILE  SIM  OUTPUT  A 
Do  you  wish  to  BROWSE  your  output?  (Y) 

Print  your  output  file?  (Y) 

n 

Do  you  wish  to  XEDIT  the  program  file?  (Y/N) 

n 

Do  you  wish  to  run  the  program  again?  (Y) 

n 

R;  T=4.63/5.82  12:11:54 
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4.  OUTPUT  1  (individual  order  policy) 


**  INPUT  DATA 


USED  DATA  FILE  NO.  =   1 
ORDER  POLICY  OPTION  =   1 
STIMULATION  LIMIT (YRS)  = 
START  LEVEL  OPTION  =   2 
PRINTOUT  OPTION  =   4 
LEAD  TIME  TYPE   =   1 
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*  ORDER  POLICY  OPTIONS 

1  :  INDIVIDUAL  ORDER  POLICY 

2  :  JOINT  ORDER  POLICY 

*  START  LEVEL  OPTIONS 

1  :  REGULAR  ( IND .  POLICY  -  E.O.Q  ;  JOINT  POLICY  -  OD.UP  QTY.) 

2  :  RANDOM  (REORDER  POINT  +  1;  ORDER-UP  POINT) 

*  PRINTOUT  OPTIONS 


PRINT  ALL  ITEM 
PRINT  ONLY  EVENT   ITEM 
PRINT  ONLY  SUMMARY 
PRINT  ONLY  FINAL  RESULT 


*  LEAD  TIME  TYPES 

1  :  COSTANT  LEADTIME 

2  :  RANDOM  LEADTIME  (UNIFORM  W/  MEAN) 


*  SERVICE  LEVEL:  P2=  99 


*  AVG.  LEAD  TIME:   0.083  YRS . (  1  MON.) 


TEM 

YR  DEMAND 

EOQ 

MUST  BUY 

CAN  BUY 

ORD.  UP 

START 

EXP.  LT 

1 

290.000 

159 

25 

112 

177 

65 

0.083 

2 

41.000 

143 

3 

35 

95 

45 

0.083 

3 

77.000 

109 

7 

48 

93 

38 

0.083 

4 

122.000 

178 

10 

77 

150 

60 

0.083 

*  INPUT  DATA  FOR  COSTS: 

ITEM  UNIT  PRICE  H.RATE  H.COST 


1 

6.90 

0.20 

1.38 

2 

1.20 

0.20 

0.24 

3 

3.90 

0.20 

0.78 

4 

2.30 

0.20 

0.46 

ORDER  COST 
GROUP  IND 
50.00  10.00 
50.00  10.00 
50.00  10.00 
50.00     10.00 


B/O  COST 
QTY.  TIME 
2.76  0.07 
0.48  0.01 
1.56  0.04 
0.92     0.02 


**  FINAL  RESULTS  (AVG.) 


DURING  50   YEARS    *** 


ITEM  #  #  DEMANDS  #  ORDERS 

1  292.24  1.84 

2  41.58  0.30 

3  77.00  0.72 

4  120.50  0.63 
TOT  AVG.;  531.32  3.54 
AVG.  NO.  OF  JOINT  REP;  0.000 


AVG.  UNIT  YEAR 
ON  HAND    B/O.S 


79.151 
69.217 

54.064 
88.760 


0.025 

0.008 
0.004 
0.014 


AVG, 
QTY. 
2.92 

0.34 
0.46 
0.84 


B/O 

TIME 
0.009 

0.006 
0.003 
0.006 


SERVICE 
LEVEL (%) 
99.00 
99.18 
99.40 
99.30 
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*  COMPUTED  AVG.  ANNUAL  COSTS  ($) 


TEM 

HOLDING 

BACKORDER 

ORDER 

TOTAL 

QTY. 

TIME 

1 

109.228 

8.059 

0.036 

110.400 

227.723 

2 

16,612 

0.163 

0.002 

18.000 

34.777 

3 

42.170 

0.718 

0.003 

43.200 

86.090 

4 

40.830 

0.773 

0.007 

40.800 

82.409 

TOTAL 

=   212.400 

430.999 

STAND. 

DEVIATION  = 

51.942 

66 


5.  OUTFUT  2   (joint  order  policy) 


**  INPUT  DATA 


USED  DATA  FILE  NO .  =   1 
ORDER  POLICY  OPTION  =   2 
STIMULATION  LIMIT (YRS)  =  50 
START  LEVEL  OPTION  =   2 
PRINTOUT  OPTION  =   4 
LEAD  TIME  TYPE   =   1 


*  ORDER  POLICY  OPTIONS 

1  :  INDIVIDUAL  ORDER  POLICY 

2  :  JOINT  ORDER  POLICY 

*  START  LEVEL  OPTIONS 

1  :  REGULAR  (IND.  POLICY  -  E.O.Q  ;  JOINT  POLICY  -  OD .  UP  QTY.) 

2  :  RANDOM  (REORDER  POINT  +  1;  ORDER-UP  POINT) 

*  PRINTOUT  OPTIONS 


PRINT  ALL  ITEM 
PRINT  ONLY  EVENT   ITEM 
PRINT  ONLY  SUMMARY 
PRINT  ONLY  FINAL  RESULT 


*  LEAD  TIME  TYPES 

1  :  COSTANT  LEADTIME 

2  :  RANDOM  LEADTIME  (UNIFORM  W/  MEAN) 


*  SERVICE  LEVEL:  P2=  99 


*  AVG.  LEAD  TIME:  0.083  YRS . (  1  MON. ) 


ITEM 

YR  DEMAND 

EOQ 

MUST  BUY 

CAN  BUY 

ORD.  UP 

START 

EXP.  LT 

1 

290.000 

159 

25 

112 

177 

71 

0.083 

2 

41.000 

143 

-1 

35 

95 

28 

0.083 

3 

77.000 

109 

5 

48 

93 

32 

0.083 

4 

122.000 

173 

7 

77 

150 

50 

0.083 

*  INPUT  DATA  FOR  COSTS: 

ITEM   UNIT  PRICE   H.RATE   H.COST 


6.90 
1.20 
3.90 
2.30 


0.20 
0.20 
0.20 
0.20 


1.38 
0.24 
0.78 
0.46 


**  FINAL  RESULTS  (AVG.) 


ORDER  COST 
GROUP       IND, 
50.00     10.00 
50.00    10.00 
50.00    10.00 
50.00     10.00 


B/O  COST 
QTY.  TIME 
2.76  0.07 
0.43  0.01 
1.56  0.04 
0.92     0.02 


DURING  50   YEARS 


ITEM  #   #  DEMANDS    #  ORDERS 


292.24 

41.58 

77.00 

120.50 


TOT  AVG.;  531.32 


1.94 
0.60 
1.18 
1.06 
4.78 


AVG.  UNIT  YEAR 
ON  HAND    B/O.S 


75.843 

56.035 
51.541 
30.055 


021 

000 
001 
003 


AVG, 
QTY. 
2.56 
0.00 
0.04 
0.30 


B/O 

TIME 
0.007 
0.000 
0.001 
0.001 


*** 

SERVICE 
LEVEL(%) 

99.12 

100.00 

99.95 

99.75 


AVG.  NO.  OF  JOINT  REP;  1.620 
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*  COMPUTED  AVG.  ANNUAL  COSTS  ($) 
ITEM   HOLDING 

1  104.663 

2  13.448 

3  40.202 

4  36.825 


BACKORDER 

ORDER 

TOTAL 

QTY. 
7.066 

TIME 

0.030 

113.400 

225.159 

0.000 

0.000 

6.000 

19.448 

0.062 

0.001 

12.800 

53.065 

0.276 

0.001 

13.600 

50.703 

TOTAL 

=   145.300 

348.375 

STAND. 

DEVIATION  = 

21.055 
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